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We present a model for the relative velocity of inertial particles in turbulent flows that 
provides new physical insight into this problem. Our general formulation shows that the 
relative velocity has contributions from two terms, referred to as the generalized accelera- 

^ ' tion and generalized shear terms, because they reduce to the well known acceleration and 

M-H , shear terms in the Saffman- Turner limit. The generalized shear term represents particles' 

C/3 ' memory of the flow velocity difference along their trajectories and depends on the inertial 

particle pair dispersion backward in time. The importance of this backward dispersion 

^ ' in determining the particle relative velocity is emphasized. We find that our model with 

^H ' a two-phase separation behavior, an early ballistic phase and a later tracer-like phase, as 

Oj' found by recent simulations for the forward (in time) dispersion of inertial particle pairs, 

gives good fits to the measured relative speeds from simulations at low Reynolds num- 
bers. In the monodisperse case with identical particles, the generalized acceleration term 

►^ , vanishes and the relative velocity is determined by the generalized shear term. At large 

0\ ' Reynolds numbers, our model gives a S'i^ " dependence of the relative velocity on the 

Stokes number St in the inertial range for both the ballistic behavior and the Richardson 

\l ' separation law. This leads to the same inertial-range scaling for the two-phase separation 

that well fits the simulation results. Our calculations for the bidisperse case show that, 
lO ' with the friction timescale of one particle fixed, the relative speed as a function of the 

^^ , other particle's friction time has a dip when the two timescales are similar. This indicates 

that similar-size particles tend to have stronger velocity correlation than different ones. 
We find that the primary contribution at the dip, i.e., for similar particles, is from the 
generalized shear term, while the generalized acceleration term is dominant for particles 
of very different sizes. Future numerical studies are motivated to check the accuracy of 

^ ' the assumptions made in our model and to investigate the backward-in-time dispersion 

of inertial particle pairs in turbulent flows. 

Key Words: Particle/fluid flows; Turbulent flows 



1. Introduction 

The dynamics of inertial particles suspended in turbulent flows is of both theoretical 
and practical importance. Its applications range from industrial processes, e.g., turbulent 
spray combustion, aerosols and raindrop formation in terrestrial clouds, and dust grain 
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dynamics in astrophysical environments such as interstellar media, protoplanetary disks, 
and planetary atmospheres. Particle collisions in a tm'bulent flow are of particular interest 
because they determine the growth of the particles by coagulation. The estimate of the 
collision rate requires the understanding of two interesting phenomena regarding inertial 
particles in turbulence, the preferential concentration and the turbulence-induced relative 
velocity. The latter is the focus of the present work. 

Our primary motivations for this study is its application to theoretical modeling of 
planetesimal formation in protoplanetory disks. Planetesimals are objects of kilometer 
size that can further grow into fully-fledged planets. The formation of these objects is 
believed t o start from the growt h of dust grains of micrometer size by collisional coagula- 



tion (e.g.. lWeidenschilling|ll980f ). Particles involved in this process have an extensive size 



range, from micrometer to kilometer, which corresponds to a range of friction timescale 
that covers all the scales (from the dissipation range, the inertial range to the outer 
scales) in the turbulence in protoplanetary disks. Therefore, a detailed understanding 
of particle collision velocities for a whole range of particle sizes induced by turbulence 
motions is crucial to investigating particle growth in these disks. Dust grain collisions in 
protoplanetary disks do not always lead to coagulation. The grains become less sticky as 
the size increases. If the relative speed is large, the collision between t wo large particles 



may cause them to fragment or to simply bounce off each other (e.g., iBlum fc Wurm I 
12008 ) . To judge the outcome of collisions between large particles thus requires an accurate 
prediction for the collision speed. 

Besides turbulence, there are other effects that can induce relative velocities between 
particles. For example, gravity gives rise to differential settling for particles of different 
sizes, which can have important contribution to the relative speed between these particles. 
In astrophysical environments, radiation pressure and the coupling with magnetic fields 
(through electric charge on the grain surface) can also play an important role. In the 
present study, we will neglect these effects, and concentrate on the turbulence-induced 
relative velocity. We aim at a physical understanding of relative velocities induced by 
turbulent motions alone, which is clearly a crucial step toward an accurate model for 
the particle collisions in the presence of both turbulence and the other complexities. 
The model presented here could be extended to take the other effects into account. 
Although our work is motivated by the problem of dust grain collisions in astrophysical 
environments , it has applications in other contexts such as droplet formation in cloud 
physics fe.g.. lSaffman fc Turneiil 19561) . 



The relative velo city of two nearby particl es induced by turbulent motions has been 



extensively studied. ISaffman fc Turned (|l956r ) considered particles with small inertia in 



the limit that the friction time, r^, of both particles is much smaller than the Kolmogorov 
timescale, r^. This limit is also expressed as St <^ 1, where the Stokes number, S't, is 
defined as the ratio of Tp to r,,. In this limit, the particle velocity at a given point can be 
approximately obtained from the (Ist-order) Taylor expansion of the particle momentum 
equation (eq. (2.1) be low) . For two partic l es wit h a separation, r, much smaller than the 
Kolmogorov scale, rj, ISaffman fc Turneij ( 19561 ) derived a result for the average radial 



relative velocity, (luirl). In the absence of gravity, there are tw o terms that contribute to 



the relative speed, a shear term and an acceleration term (see lAvala et g^.l 12008.) . 
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\ 1/2 
{\Wr\) = \j^ [j^-r' + (Tp2 - Tplfa^J (1.1) 

where e and v are, respectively, the average dissipation rate and the kinematic viscosity 
of the turbulent flow, Tpi and Tp2 are the friction timescales of the two particles, and a 
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is the rms acceleration, i.e., a^ = ( (■§^) ). The factor \/2/7r is from the conversion of 

the radial relative velocity variance, {w'^), to (Iwrl) assuming a Gaussian distribution for 
Wr- Note that the shear term is from the longitudinal structure function, Sii{r)^ of the 
flow, which is given by j^j^r^ for r ^ 77. In the SafFman- Turner (S-T) limit, the particle 
velocity follows the local flow velocity very closely, thus the velocity of the two particles 
at a distance r < 77 is highly correlated because of the strong flow velocity correlation 
across a small distance. 

The opposite limit is that of heavy particles with friction timescales much larger than 
the Lagrangian correlation timescale, T^. In this limit, the velocities of two particles 
coming together are essentially uncorrelated. This is because particles with Tp ^ Tl 
have long-time memory, and their current velocities have substantial contributions from 
the flow velocities on their trajectories in the past. These flow velocities at early times 
are likely to be uncorrelated because the particles were far away from each other. The 
relative velocity in this limit is thus determined by the sum of the velocity variances of 
the two particles, 

-,1/2 

(1.2) 



{\Wr\) = 



(^;'(l 



)\ +(x;'(2) 



where w'^^^ an d v'^"^^ denote the p article rms velocities. In the limit Tp ^ Tj,, they are 



given by (e.g.. lAbrahamsonlll975l ). 

(y'i^)Y^u'^lL. L'(^)Y^u''I± (1.3) 

V / Tpi V / Tp2 

where u' is the rms of the flow velocity fluctuations. In the derivation of eq. (1.3), the 
temporal correlation of the flow velocity on a particle's trajectory is approximated by the 
Lagrangian correlation for tracer particles. Discussions on the validity of this assumption 
will be given in §2.2. 

The problem of the relative velocity of inertial particles in these two extreme lim- 
its is physically clear, and the results given above are easy to understand and are ex- 
pected to be generally robust. On the other hand, for particles with intermediate inertia, 
T,j ^ Tp <Tl, the problem is more complicated and is less well understood. The velocities 
of two nearby particles with intermediate Tp are partially correlated and the degree of 
correlation, intermediate between the two limits, is not easy to evaluate. We will point 
out that a very important factor in determining this correlation is the distance between 
the trajectories of the two particles as a function of time before they come close to each 
other. To our knowledge, this point has not been clearly recognized or explicitly empha- 
sized in the literature. We will show how the separation of two nearby inertial particles 
backward in time affects the relative velocity between particles with intermediate inertia, 

'^ri r^ "^p r^ -^ L- 

A successful theory for particles of all sizes needs to correctly predict and explain the 
behavior of the relative velocity between particles with intermediate friction time, as 
well as recover the results in the two extreme limits. An example of particular theoretical 
interest is that of identical particles, referred to as the monodisperse case. In the S-T 
limit {Tp <^ T,i), the acceleration term in eq. (1.1) vanishes for the monodisperse case 
and the relative velocity does not depend on the friction time. It is constant at a given 
distance, r, and increases linearly with r. In the opposite limit with Tp ^ Tl, the relative 

velocity decreases with the friction time as Tp or St~^/'^, according to eqs. (1.2) and 
(1.3). The question of how {\wr\) scales with Tp for Tp in a range corresponding to the 
inertial range of the turbulent flow, and how it connects with the two extreme limits has 
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not been systematically studied or fully understood. This is one of the primary goals of 
the present paper. We find that the dispersion of particles backward in time is crucial to 
answer this question. 

The existing models have very d ifferent predictions for the relative speed in the inertial 
range for the mo nodisper se case ( Volk et al. 1980; Williams & C randll983t lYuulll984 : 
Kruis fc KusterdlT997: Zaichik et al 2003, 2006: Avala et al. 2008) . A detailed discussion 
of the qualitative differences between these models and their problems will be given in §4. 
An important reason for the problems in most of the previous models is that they did not 
clearly recognize or carefully consider the eff ect of the part i cle pa i r sepa ration backward 
in time (except for the differential model bv lZaichik et al\ ( 20031 120061) to be discussed 
below, which we think has the particle backward separation implicitly included). The 
role of this backward separation will be revealed along the formulation of our model. 

In the previous studies that cover a whole range of Stokes numbers, the differential 
mode l by Zaichik and collaborators ([Zaichik fc Alipchenkovl l2Q03t IZaichik et a/.l 12003 , 
20061 ) is perhaps the most complete one, as it examines the effect of preferential cluster- 
ing and the relative speed simultaneously. We will refer to this model as Zaichik et aVs 
model or the model by Zaichik et al. Assuming Gaussian statistics for the flow velocity, 
the model first sets up an equation for the joint probability distribution function (pdf ) of 
the particle separation and the relative velocity. Deriving the first 3 moment equations of 
the pdf equation and closing these moment equations by a quasi-normal approximation, 
Zaichik et al. were able to obtain a set of partial differential equations for the particle 
density correlation function (usually referred to as the radial distribution function) and 
the particle velocity structure functions. The solution of the differential equations repro- 
duces the two extreme limits discussed above and predicts that the relative velocity of 
identical particles with intermediate inertia is proportional to Tp or, equivalently, St^''^ . 
The validity of this prediction remains to be confirmed by high-resolution simulations. 
Despite the elegant mathematical formulation and good agreement with results of direct 
numerical simulations of turbulence with low Reynolds numbers, the model lacks physi- 
cal transparency in its approximations, especially the quasi-normal assumption. One of 
the goals of our model is to elucidate the physics behind the inertial range scaling of the 
relative velocity. 



Falkovich et al. I (|2002f ) found that an e ffect, named the sling effect , has a significant 
contr i bution to the relative spee d (see also IWilkinson fc Mehlig Il2005t IWilkinson et al. 



2006t iFalkovich fc Pumir 1 120071 ) . The physical picture of the effect is that, at regions 



with large negative velocity gradients, fast er moving particles can catch up the slower 
ones f rom behin d ( Falkovich et al. 1 12002|) , leading to trajectory cros sing of the par- 



ticles (jBec et al\ (|2005l ), also see Fig. 1 in IFalkovich fc Pumir I (l2007|) for an illustra- 



tion). This results in a lar ger relative speed than the S-T prediction for small particles. 
Falkovich fc Pumir I (J2007I ) showed that the effect starts to be important for St <; 0.2 and 
gives a relative speed several time larger than eq (1.1) for St between 0.2 and 1. We will 
point out a common element shared by the sling effect and our model: the contribution 
to the relative speed from the particles' memory of the flow velocity difference in the 
past. 

The paper is organized as follows. We present the formulation of our model in §2 (a 
general formulation in §2.1, and basic assumptions and approximations in §2.2 and §2.3). 
The results for identical particles (monodisperse) and different particles (bidisperse) are 
given in §3.1 and §3.2, respectively. In §4, we discuss previous models and compare them 
with our model. Conclusions are given in §5. 
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2. The Model 

2.1. The General Formulation 

The velocity, v{t), of a particle with friction time, Tp, can be obtained by integrating the 
momentum equation, 

dv _ u{X{t),t) -V 



dt Tp 



(2.1) 



where u{x,t) denotes the flow velocity field and X{t) is the position of the particle at 
time t. Clearly, u {X{t),t) is the flow velocity at the positions of the particle along its 
trajectory (we will refer to it also as the flow velocity "seen" by the particle). The particle 
trajectory is given by, 

X{t) =Xo+ [ v{t')dt' (2.2) 

Jta 

where Xq is the particle position at a given time ig- 
Equation (2.1) has a formal solution. 



vit) = t)o exp ( ] ^ I u {X{t), t) exp ( | dr 

V ^P / ^P J to \ '^P J 



(2.3) 



where Vq is the particle velocity at io- 

We are interested in deriving the relative velocity between two particles at a distance 
r at a given time t. We label the two particles by superscripts "(1)" and "(2)". For 
example, their velocities at t are denoted as v'^-^^t) and v'^^^{t), respectively. To evaluate 
the average relative speed, we will calculate the velocity structure tensor, Spij, of the 
two particles, 

where (• • •) denotes the ensemble average. The particle velocities can be solved by inte- 
grating equation (2.3) and the trajectories of the two particles are subject to a constraint, 

X^^'>{t)-X^^'>{t)=r. (2.5) 

which means that two particles happen to be separated by r at t. We will particularly 
consider small values of r (below rj) because we are interested in the collision speed, which 
is essentially the relative velocity of the two particles over a distance of the particle size. 
From the structure tensor Spij, we will obtain the longitudinal structure function Spu, 
which, by definition, is the radial relative velocity variance, {w'^). Although only small 
r will be considered in the paper, our model can predict the structure function at all 
separations. By a comparison with their results, ou r model may provi de an explanation 



for the inertial particle structure functions found in lBec et al\ (|2009bl ). 
The particle structure tensor can be written as, 

^P. = {v^^v^') - {vi'^'^) - {vM') + (vPvf) . (2.6) 

Note that the cross terms correspond to the particle velocity correlations discussed in §1, 
where it was argued that a careful treatment of such correlations is essential for modeling 
the relative velocity of particles with intermediate inertia. 

To calculate the structure tensor, we insert eq. (2.3) for the particle velocities into eq. 
(2.6). For simplicity in notations, we will set the time when the particle relative speed is 
measured (i.e., t in eq. (2.4)) to be zero, and assume it is far from the initial time (since 
we are interested in the relative velocity for a statistically stationary state). This allows 
us to set to in eq. (2.3) to — cxj. 
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We analyze the four terms in eq. (2.6) one by one. The first term on the r.h.s. corre- 
sponds to the velocity variance of particle (1). For this term, only the velocity of particle 
(1) is involved and we have, 



where 



vf^4^\ 



,(i)m- 



— I — lu^\T)uf\r'))e^^ ( _ ) exp ( — 



Tpl 



Tpl 



Tpl 



T 
Tpl 



(2.7) 



'(i) 



'(1)/ 



X^ '{t),t\ denotes the flow velocity on the trajectory of particle (1). 



The exponential factors here represent the memory loss of the particles. The integral 
limits in eq. (2.7) (see also eq. (2.9)) suggest that it is the flow velocity the particles 
saw in the past that is relevant in determining the particles' velocities at the current 
time. The relative position of the two particles back in time will be shown to play an 
important role in the prediction of their relative velocity. We will call ( u| {t)u)- (t') 

in the integral the trajectory correlation tensor and denote it as -B-j/,-, i.e.. 



B^,ir.r') = iv^^(r)^^ir')) 



where the subscript "T" stands for "trajectory" . 

The result for the 4th term on the rhs of eq. (2.6) is similar. One only needs to replace 



j(i) 



,(2) 



,(2)/ 



,(2), 



{t)u, (t')). If the two particles are 



Tpl in eq. (2.7) by Tp2, and B!^-'j by B^/^- 

identical and have the same friction time, (vl'^'vY') is equal to ivY'vY')- These two 

terms correspond to the velocity variance of each particle and will be called the velocity 
variance terms. 

The exact form of B^ij as a function of the friction time is not available. In the 
limit of vanishing Tp (i.e., passive tracers), this correlation tensor would approach the 
Lagrangian correlati on tensor, -BLij, of the flow, which has been extensively studied (e.g., 
Yeung fc Popelll989). A common approxiination i s to set Btja equal to B^,i^ for particles 



with any Tp (e.g., IZaichik fc Alipchenkovl l2003t IZaichik et al\ l2003l l2006t lAvala et al 



20081) . Physically, it corresponds to the assumption that the trajectory of any inertial 



particle is not far away from that of a tracer particle starting from the same initial 
condition. We will adopt this assumption in our calculations and its validity will be 
discussed in §2.2. 

The cross correlation terms in eq. (2.6) can be evaluated with the same approach. The 
second term on the r.h.s. is given by, 

'•0 W^ /-o „ , , , 

exp ( — ) . (2.9) 



(1) (2) 



dT 

Tpl 



dr^ 

Tp2 



(.«(r)uf(r'))exp(-^) 



T 
Tp2 



(2) (1) 



in eq. (2.6) is similar to eq. (2.9). The sum of these 



The result for the term 

two tensors can be written as, 

(u^\T)ufiT')) + (u'^\r)uf\T')) = 4V,.(r,r') +4'i(^,^') " ^T., (r; r, r') (2.10) 
where the tensor Stij is defined as 



STij{r;T,T') 



,(1) 



(r) 



,(2) 



(r) 



,(1) 



(r') 



,(2) 



(r') 



(2.11) 



Clearly, Stij is the correlation of the flow velocity difference at the positions of the parti- 
cles on their trajectories at times r and r'. We have explicitly indicated the dependence 
of the tensor on the particle separation at time zero. The ensemble average on the r.h.s. 
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includes an average over the probability distribution of the flow velocity difference at 
time zero. By analogy with -Bxij, we will call STij the trajectory structure tensor. 

Since STij has not been directly studied, we will give an approximate estimate for it 
in §2.3. For example, we will relate the flow velocity difference, u^^\t) — u^^'>{t), along 
the trajectories by the two particles, to the separation of the two particles at t, assuming 
the velocity difference scaling in the Eulcrian frame applies to the velocity difference on 
the particles' trajectories. The uncertainty in this assumption will be discussed in §2.3. 
We will denote the particle separation at a given time t as p(t), which, in our notation, 
is given hy X^ ' (t) — X^ ' (t). Note that p is a stochastic vector because of the particle 
dispersion by turbulent motions. 

Combining eqs. (2.6), (2.7), (2.8), (2.9) and (2.f0), we finally arrive at the formula for 
the velocity structure tensor of two particles separated by r. 



Sp^J{r)=A^J+V,J, (2.f2) 



where 



J, /-O 



J-oc>Tp2 J-ooTp2 ^ \Tp2j \Tp2 



T 



and 

'^v^ I —f — 5T»j(r;T,T')exp( — )exp( — ). (2.f4) 

J-ooTpl J-ooTp2 V'^pl/ \Tp2/ 

In the trajectory structure tensor in "Dij, the dependence on r is from the requirement 
that particle separation p is equal to r at time zero, i.e., 

p(0)=r. (2.f5) 

The result, eq. (2.12), is written in such a way that the first term Aij only depends 
on the 1-particle trajectory correlation tensor and the second term Vij only on the 2- 
particle trajectory structure tensor. There are also physical reasons to split Spij into 
these two terms. First, the Aij term vanishes for identical particles with Tpi = Tp2, and 
only Vij contributes to the relative speed in the monodisperse case. On the other hand, 
for particles of very different sizes, Aij dominates the contribution to the relative speed 
(see §3.2). Second, in the S-T limit, Aij and Vij reduce to the acceleration term and the 
shear term in eq. (1.1), respectively. Therefore, our formulation can be regarded as one 
that extends eq. (1.1) from the low-inertia limit to the whole range of particle sizes. We 
will refer to Aij and Vij as the generalized acceleration term and the generalized shear 
term, respectively. 

It is straightforward to see that Vij reduces to the shear contribution in the the S-T 
limit. As Tp — )■ 0, we have l/Tpexp(T/Tp) — > S(t), and thus Vij approaches the flow 
structure tensor Sij{r). As pointed out in §1, it is exactly this flow structure tensor that 
is responsible for the shear contribution in eq. (1.1). 

In the S-T limit, Aij can be evaluated as follows. For Tp — > 0, particle trajectories 
are close to those of tracers, so i?Tij — Bi^ij- For small time lag At = t — t' (only 
small time lag is of interest here because of the exponential cutoffs in the integrand), 
Bi^ij{/S.T) ~ (m'^ — ^a^AT^)(5ij where u' and a are, respectively, the rms flow velocity and 

the rms acceleration. Using this approximation for both B\.-- and Bi^^- in eq. (2.13), we 
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T-p2?h 



which is exactly the same as the acceleration term in eq. 



find Aij — a^(Tpi 

(l-l)- 

In §2.2 and §2.3.3 we will show that, with our modeling of Aij 

(2.13), and (2.14) also recover the large particle limit as well. 

Our formulation reflects the trajectory-crossing effect mentioned in §1. When setting r 

to zero (or more exactly the particle size) , the formulation is for two p articles whose tra- 



and 2?ij, eqs. (2.12), 



jecto ries cross at time zero. Different from the model for the sling effect (jFalkovich et al. 



20021 ) , our model does not specify the physical mechanism how and when the trajectories 
of two particles cross. Instead the process leading to trajectory-crossing is considered in- 
directly from the backward separation behavior of the two particles. From the following 
perspective, one may see a common feature shared by the sling effect and our model. The 
sling effect could be interpreted as a mechanism that contributes to the particle separa- 
tion backward in time, and it gives a larger relative speed by increasing the contribution 
from the particles' memory of the flow velocity difference in the past. The latter is the 
point of our model. Therefore, we think that the sling effect can be accounted for in our 
model if the backward separation to be used in the model includes its contribution. 

2.2. Modeling Aij 

The formulation has been general so far. To proceed, we make assumptions for the tra- 
jectory correlation and structure tensors in eqs. (2.13) and (2.14). In this subsection, we 
evaluate the generalized acceleration term Aij . 

We use the usual assumption for BTij that the flow velocit y viewed by a particle 
on it s trajectory is the sa,me a s that by a tracer particle (e.g.. IZaichik fc Alipchenkovl 



2003t IZaichik et a/.|[2003l l2006t lA^ala et a;.ll2008h . i.e., Bt,j = Bi^.j. In a statistically 



stationary and isotropic flow, the Lagrangian correlation tensor only depends on the 
time lag, and can be written as Bi^ij{/S.T) — u''^5ij^{/S.T) where $ is the normalized 
temporal correlation function. We adopt the bi-exponential form for $ (see ISawford 
199lUZaichik et adl2003ll2006l ). 



$(AT;rT,TL) = 



1 



2Vl - 2z2 



l^l + Vl - 2z2) exp 



2|A7 



(l - Vl - 2z2) exp 



(1 + VI - 2z2) Tl 
2|Ar| ^ 



[l - Vr^2i^) Tl 



(2.16) 



where Tl {— / <I'(Ar)(iAT) is the Lagrangian correlation timescale, and z — tt/Tl is 
the ratio of the Taylor micro timescale, tt, to T^. In addition to the time lag Ar, we 
have written $ also as a function of tt and Tl in eq. (2.16) for later convenience (see 
eq. (2.23)). 

We obtain Tl using DNS results for the ratio of Tl to the large-eddy turnover time, 
Te- Te is defined as the longitudinal integral scale, Li, divided by the rnis velocity, u' . 
The length scale Li can be calculated from the relation, e = Du'^ /Li, using simulation 
results for the dimensionless coefficient D. Defining a length scale L as L — u'^ /e, we 
have Li — DL. The large-eddy ti mescale is then given by Te = Du''^ /e or Te = DT^, with 



Te defined as T 
is 



x'2/e = L/i 



Yeung et al\ (|2006l ) found that the ratio of Tl to Te 
0.75, and is essentially independent of the Taylor Reynolds number. Rex (although 
this ratio may depend on the forcing and thus could be flow- dependent). Therefore we 
have Tl = 0.75Du'^/e or Tl ^ 0.75DTe. 

Numerical simulations have shown that the coefficient, D, is Reynolds number depen- 
dent for Re\ ^ 10 0, but approaches a constant ~ 0.4 for Re\ larger than several hundred 
(e.g.. lYeung et al...200(k Is hihara et qL,2009,) . Therefore in the limit of large Re\, we ex- 
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pect Te = 0.3u'^/e. To also account for the Rex dependence of D at small j^cai we use 
P = .4(l + 30/i?eA), which is obtained from fitting the numerical results in lYeung et all 
(|2006[ ). We then have, 

Tl = 0.3(1 + SO/Rexy^/e (2.17) 

which is very cl ose to that ad o pted by IZaichik et all (|2003l 120061 ) and the empirical 
formula given in ISawford et all (|2008l ). The definition of Rex gives m'^ = Rex/V^u'^, 
which will be used for normalization in our calculations. 

The Taylor micro timescale is defined as tt = (2u'^/a^)^/^. The asymptotic behav- 
ior of the normalized acceleration variance, ao — d} l{^l'^v^^l'^\ at large Rex has not 
been resolved by current simulations. Although g p is predicted to be constant at larg e 
Rex limit by the Kolmogorov 41 theory (see, e.g.. lVoth et a^.lll998tlZaichik et aZ.ll2003l ). 
the intcrmitt ency correctio ns to the K41 theory may give it a power-law dependence 
on Rex fe.g.. lBorgaslll993l ). For example, assuming that the temporal statistics of the 
dissipation rate along Lagrangian trajectories are the same as its spatial statistics in 
the Eulerian f rame (which follows from the ergodic hy pothesis and incompre ssibility; see 



BorgasI ( 19931 ) ) , and using the intermittency theory by IShe fc Leveaud ( 19941) for the dis 



PC 

1 



sipatio n rate statistics, we find that aq oc i?e^'^"^'^ (a similar result was obtained bylBorgas 



(|1993I ) using the log- normal intermittency model) . This result is in impressive agreement 
with one of two formulas th at well fit the results from simulations with resolution up to 



2048^ in lYeung ei oil (|2006[) 



ao = 1.9i?e5-i3^(l 



85/i?e 



which goes like R^^^^^ at large Rex- Therefore, one may expect that oq 



(2.18) 

1.9i?e5"5 

for asymptotically large Rex- However, the confi rmation of this as ymptotic behavior 
needs future simulations wit h higher resolut ions ( Yeung et aLl 120061 ) or more accurate 
experimental measurements ( Voth et aLlll998l ). We will use eq. (2.18) in our calculations. 
Eq. (2.16) approaches exp(— |Ar|/T£,) for z <C 1 and |Ar| '^ tt- The biexponential 
form is expected to be better than the single exponential form, exp(— |At|/Tl), because 
the Lagrangian correlation is believed to be smooth at small time lag, |At| — > 0. It can 



1 - AtVt. 



^ for 



be easily shown that $ 
requirement. 

We insert eq. (2.16) into eq. (2.13) to calculate Aij 
integration gives. 



Ar| ^ Tt, and thus satisfies the smoothness 



A lengthy but straightforward 



y\i 



u"5, 



{^2-^1^ [Vlin2 + {^l+^2)- 



{n^ + n^) \vii + Q.i + — 



^2 



^1 + ^ 



(2.19) 



where Qi and ^2 are defined as fli = Tpi/T^ and U,2 = Tp2/TL. Equation (2.19) correctly 
reproduces the acceleration term in the S-T limit. When fii, ^2 ^ 2:^/2, i.e., Tp^, Tp2 <C 



T|./(2Ti) (which is '-^ t,, from the K41 phenomenology) , we have Ai 



f]l)V(zV2) = {Tp2 



Tpi^a'^S,^ 



H,j{n2 



where the definitions of fi, z and tt have been used in 



the last step. As expected, for identical particles, the Aij term is zero. 

Results closely rel ated or directly comparab le to eq. (2.19) have been derived from 
several other models. IWilliams fc Grand (1983) considered the relative velocities in two 
limits with 17 <C 1 and fi ^ 1, and then gave a "universal" solution by interpolating 
those two limits. For the small particle limit, they argued that the particle separation 
(back in time) can be neglected in the calculation of the particle velocity correlation as 
long as one of the two particles is very small. With this assumption, they find that the 
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relative velocity variance, (w^), is given by (their eq. 19), 

This result corresponds to our result for Aij because neglecting the backward separation 
is essentially the same as neglecting the 2?^ term. The latter can be justified if one of the 
two particles, say particle (1), has a very small friction time, i.e., Tpi -^ 0. In that case, 
eq. (2.14) is approxim ately given by f ^^ STij{r;0,T') exp{T' /Tp2)dT' /Tp2. For r ^- as 



considered by William s fc Grand (|l983l ). we have STij{r,0,T') — >■ and thus Vij — >• 0. 

It is obvious that, if z is set to zero, eq. (2.19) for Aij reduces to eq. (2.20). This 
is expected because, wit h z = 0, $ take s the same single-exponential form used in the 
derivation of eq. (2.20) bv I Williams fc Cr ane (198^). Without the z terms, eg . (2.20) does 
not reproduce the acceleration term in th e S-T limit. iKruis fc KustersI ( 19971 ) generalized 



the model by IWilliams fc Grand ( 1983 ) and used a temporal velocity spectrum that 



accounts for the acceleration field in the flow (corresponding to a correlation function 
similar to our eq. (2.16) ) and obtained a formula that gives the acceleration term in eq. 
(1.1) in the S-T limit. 

In §3.2, we find that for two particles of very different size the contribution to the 
relative speed from Aij dominates over that from Vij, and thus eqs. (2.19) or (2.20) can 
be used to estimate the relative velocities between very different particles (although the 
justification above for neglecting Vij is only for the case with at least one tiny particle). 
However, for similar particles, Vij is the dominant term, and thus the result for J7 <C 1 
by Williams and Grane is not valid for similar-size particles. In the limit of particles 
with 17 ^ 1, Will iams fc Grand ( 1983 ) considered particle separation. We will give more 



comments on their model in 52. 2. 3 and §4. 



Yuul (|1984l ) derived a formula for the relative velocity with a shear term and an acceler- 
ation term. The derivation included an added mass term, bdtu, in the particle momentum 
equation (eq. (2.1)). The coefficient b, given by 3pf/{2pp + pf), is small for solid particles 
in a gaseous flow and the added mass effe ct (K r uis an d Kuster 1997) is not important. 
In that case the acceleration term given bv I Yuul (|l984l ) is exactly the same as eq. (2.20). 
This means that the acceleration term by Yuu is generally in agreement with our eq. 
(2.19) for Tp 3> t^/Tl ~ Tjj. His shear term thus corresponds to our Vij term, which 
is the dominant term for similar particles. The problem in Yuu's shear term is that its 
derivation did not keep track of particle distance in the past, and thus the resulting 
shear term does not account for particle's memory of the flow velocity difference. This 
would underestimate the relative velocity between similar particles since the flow velocity 
difference was larger at earlier times when the particle separation was larger. 

With the Lagrangian correlation function, eq. (2.16), the flrst and third terms in eq. 

(2.13), representing the velocity variances of particle (1) and (2), (^v'^^^) and (u'^^-*) , 
are given by. 



"'-")" = »" "-r.^:. . ...<i (.■'")' ^ «" ''-:'''.. ■ (2.21) 



2 

ni + nj + z^/2' ' V J n2 + ni + z^/2' 



For r^i, 1^2 "^ 1, the particle rms velocity is close to the flow rms velocity u' . As expected, 
in the limit ili, ^2 ^ 1, the sum of the two terms reproduces eqs. (1.2) and (1.3) for the 
relative velocity of large particles. This means that, in order to build a model that gives 
correct prediction in the limit of very large Tp, one needs to guarantee that the sum of 
the other two terms in the model, i.e., the 2nd term in eq. (2.13) and Vij (both from 
the cross correlation of particle velocities), approaches zero as Tp —>■ oo. More speciflcally. 
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the sum of those two terms has to approach zero faster than l/r^ so that they do not 

dominate over (w'^^^) and (u''^') given by eq. (2.21). We will show that the model 
presented in the next subsection for the Dij term does satisfy this constraint, and hence 
eqs. (1.2) and (1.3) are correctly reproduced in our model. 

Finally, we point out that the assumption we adopted that the temporal correlation 
of the flow velocity on an inertial particle's trajectory can be approximated by the La- 
grangian correlation function may be invalid for large Tp. The trajectory of a large particle 
can be very different from that of a tracer particle. For example, if Tp <; T^ or Te, the 
particle may not move significantly as the flow "sweeps" by. Thus the correlation of the 
flow velocity along the trajectory of such a heavy particle may be better approximated 
by the Eulerian temporal correlation. This means that the flow velocity correlation on a 
particle's trajectory could make a transition from Lagrangian-like to Eulerian-like as Tp 
increases. In that case, replacing the Lagrangian correlation timescale Tl in eq. (2.21) by 
the Eulerian correlation timescale, Tem, would give a b etter estimate for the rms velocity 



of large particles. Tl in eq. (1.2) from lAbrahamsonl ()l975[ ) shoul d also be replaced by 
Teu- If the Eulerian co rrelation timescale is larger than Tl (e.g., lYeung fc PopdllQSil 
Kaneda fc Gotoh 1991). our model would underestimate the relative speed by a factor 



Kaneda fc Gotoh 1119911 ). our model would underestimate the relative speed by a factor 
of {Teu/TlY^'^ in the limit of Tp ^ T^. A numerical study of the flow velocity correlation 
on particles' trajectory as a function of Tp would be useful to improve our model. 

2.3. Modeling Vij 

In order to evaluate the generalized shear term, I?^, we need to model the trajectory 
structure tensor, S'tij, which has not been directly measured. The tensor can be formally 
written as 

St.,^ J J {Su'f\T)Suf{T')\p,p')p{p,p';r,T,T')dpdp' (2.22) 

where P is the joint probability distribution of the particle separations, p and p', at 
r and r', respectively, and we have used (5m^^ (r) (= Ui(X'-'-^(r)) — Ui(X'^)(T))) and 

5uf \t') to denote the flow velocity difference "seen" by the two particles. The ensemble 
average term in the integrand is the velocity difference correlation conditioned on the 
particle separations. A series of assumptions need to be made for the estimate of S'tzj 
since both the conditional correlation and the joint probability of the particle separations 
are unknown. 

The conditional correlation depends on both the separations, p and p' , and the two 
times, T and t'. It may also have a direct dependence on r, in addition to that through 
p and p' . This possible dependence is neglected here. We assume that the amplitude 
of the velocity difference 5uf (r) on particles' trajectories can be approximated by the 
Eulerian velocity difference, Sui{p,T), across th e sepa. ration p. This assumption is sim- 
ilar to the ind ependence hypothesis bv ICorrsinI ( 19591) (note the simi l arity b etween eq. 



(5) in ICorrsin (1959) and our eq. (2.22); see also IShhen fc Corrsin"! (|l974l) ). Assum- 
ing that the displacement of a fluid particle (f rom its i nitial position) is statistically 
independent from the particle's current velocity, [Corrsid ( 19591 ) gave a relation between 



the Lagrangian correlation function and the Eulerian correlation. The hypothesis by 
Corrsin essentially neglects a constraint between the particle displacement and the ve- 
locity along the Lagrangian trajectory (i.e., the latter equals the time derivative of the 
former). In the case of inertial particles, the constraint neglected in our approximation 
is d'^p/dt^ = {Su^P^ — dp/dt)/Tp. The existence of this constraint could give rise to a 
statistical correlation between Su'^'p^ and p. For example, a larger particle separation at 
a given time suggests a larger (on average) velocity difference in the past, and perhaps 
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a higher possibihty of a large velocity difference right at that time. Therefore, relative 
to the Eulerian velocity difference across a fixed separation, the velocity difference along 
the particles' trajectories may have a stronger dependence on the particle separation. 
Setting 5uf' to the Eulerian velocity difference, 5ui{p), could thus underestimate the 
dependence of Suf on p. The uncertainty and reliability of our assumption here are 
subject to tests by future numerical experiments. 

From the assumption above, we have {du^' {t)Suj {t')\p, p') ~ {5ui{p,T)5uj{p' ,t')). 
The correlation of the Eulerian velocity differences at two times is also unknown and 
further approximations are needed. We first assume that it can be written as a product 
of a separation dependence term and a time-lag dependence term. The separation de- 
pendence term is then assumed to take the form of the Eulerian structure function of the 
flow, Sij. Note that, the conditional correlation actually depends on two separations, p 
and p' , but for simplicity we will approximate it by Sij{R), at a single separation, R, 
characteristic of the particle distances between r and r'. The choice for the separation 
ii as a function of p and p' will be discussed later. As a function of the two stochastic 
vectors, R is also stochastic. 

The time-lag dependence accounts for the temporal correlation between the flow ve- 
locity differences "seen" by the two particles. This correlation depends on the persistence 
of the "structure" in question, which is a function of the "structure size" . Here the size is 
essentially the distance between the two particles. Associated with each structure, there 
is a correlation timescale, Tr. To estimate Tr, one needs to pick up a size to characterize 
the structure that corresponds to particle distances at time between r and r'. We will 
take the size to be the same as the particle separation R to characterize the separation 
dependence. We point out that there is no physical motivation for this particular choice, 
and that it could be better to choose the distance, Rm, at the earlier one of the two 
times T and t', min(r, r'), assuming the persistency of a structure is determined by its 
initial size. The difference in the results from the two choices will be discussed in §3. 
The time lag dependence is assumed to take the same function form as eq. (2.16) for 
the Lagrangian temporal correlation in §2.2. Namely, we set the time-lag dependence to 
be $(t — t';ttr{R),Tr{R)), where the Lagrangian timescale, T^, in eq. (2.16) has been 
replaced by T]i{R) and tt by ttr{R)- The tii nescale ttrJR) is an analo gue of the Taylor 
micro timescale for a structure of size R fsee lZaichik et aLll2003ll2006l ). 



The conditional correlation is now approximated by Sij{R)^{T — r'; ttr{R), Tr{R)^ . 
In order to estimate the trajectory structure function, in principle one needs to integrate 
the conditional correlation over the distribution P{p,p') (eq. 2.22), or equivalently over 
the distribution, P{R)^ of R (in our approximation the conditional correlation depends 
on the separations p and p only through R). We take a simple approach here. We set the 
integral to be equal to the conditional correlation at a distance corresponding to the rms 
of R. This means that the particle distance at a given time is assumed to be single- valued, 
i.e., the distribution P{R) is a delta function at the rms of (R). A rough justification of 
this approximation will be given in §2.3.3. The direction of R probably has a random 
distribution due to the turbulent dispersion, and we will average the structure tensor 
Sij{R), over the direction distribution of R. 

For simplicity in notations, hereafter we will use R to denote the rms length of R, i.e., 
R^ = (R^) (note that this is different from the conventional notation that R denotes 
the length of a vector R, i.e., R^ = R^ , without ensemble averaging). Similarly, p(r) 
will denote the rms length of p{t). We refer to R and p as "rms distance" or simply 
"distance" , while using "separation" for the corresponding stochastic vectors. 
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The approximations for the trajectory structure tensor are now complete, 

ST^J{r;r,T') c^ {S,,{R))^^^^t ^ r':,rTRiR),TniR)) (2.23) 

where the ensemble average for 5*.^ (i?) is over the direction distribution of R. The angular 
average will be carried out in §2.3.3. 

Finally we need to specify the rms distance i? as a function of p{t) and p{t'). If the 
flow velocity difference scales with the distance as a power law, which is probably the 
case in well-developed, homogeneous and isotropic turbulence, a good choice would be 

R{r,r')^{p{T)p{T')f\ (2.24) 



We no te that, in their assumption for the Lagrangian structure tensor. I Zaichik fc Alipchenkov 
(|2003l ) apparently set R to be the particle separation at the earlier time of r and r'. We 
argue that eq. (2.24) is probably a better assumption because Snj is expected to be zero 
if either p(r) or p(t') is zero. 

The assumption, eq. (2.24), is expected to be valid when p and p' are in the same 
length scale subrange since in that case the scaling of velocity difference across the two 
separations follows the same power-law (see, §2.3.1). On the other hand, if p and p' are 
in different subranges, the structure function across R defined by eq. (2.24) may not 
correctly represent the product of the velocity difference amplitudes across p and p'. 
Fortunately, we find that this does not significantly affect our prediction of the relative 
speed based on a consideration of the $ term in the trajectory structure tensor (the 
temporal correlation of the velocity differences) . The $ term gives an exponential cutoff 
when the time lag is large, i.e., for very different r and r'. Note that, if p and p' are 
in different subranges, in general r and r' would also be very different (although it is 
possible that the particle pair experiences a large separation change during a short time 
interval, these extreme events should be rare, and not affect the low order statistics, i.e., 
the second order particle structure function, we study here). Therefore, the exponential 
cutoff from the $ term would suppress the contribution from very different p and p' to 
the integral for T^ij. In other words, the temporal decorrelation of the velocity difference 
over large time lags suggests that the main contribution to Vij is probably from similar 
p and p', where eq. (2.24) is valid. 

2.3.1. The Flow Structure Tensor and the Timescales 

In homogeneous and is otropic tur bulence, the Eulerian structure tensor at a separation 
I can be written as (e.g.. lMonin fc Yaglom 1975), 

S^j{l) = ^„„(0% + {Sn{l) - 5„„(Z))^ (2.25) 

where Snn and Su are, respectively, transverse and longitudinal structure functions. For 
an incompressible velocity field, they are related by, 

Snn{l)=Su{l)+'^^^. (2.26) 

The longitudinal structure function, Su, in different ranges of length scales is given as 
follows. 

As mentioned in §1, in the viscous subrange, Su is given by, 

Su{l) = -^l^ for l<,v. (2.27) 

In the inertial subrange, we have, 

Su{l) ^ C{eiy^\ forry<SZ^Li (2.28) 
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where the coefficient C for the v elocity scaling in the inertial range is beheved to b e 
universal and will be set to C = 2 ( Monin k, Yaglomlll975tlZaichik fc Alipchenkovll2003l ). 
For / larger than the integral scale, Sa is constant, 



5»(0 = 2m'2, for l^Li. 



(2.29) 



The characteristic scale at which Su switches from the viscous-range scaling to the 
inertial-range scaling can be obtained by equating eqs. (2.27) and (2.28). This gives a 
transition scale of (ISC)'^ /'^??, which is about I Sr? for C = 2. It is consistent with the 
simulation results given in llshihara et al\ ( 20091) . where the switch occurs at about a few 
tens of Kolmogorov scale, between 77 and the Taylor micro scale A. The scaling changes 
from eq. (2.28) to eq. (2.29) at I ~ (2/C)^/ ^L where L = M ^^/e, as defined earlier. This 
is also in general agreement with Fig. 7b in llshihara et al\ ( 20091 ). which shows that Su 
becomes constant at ~ 0.5 — iLi (Li ~ Q AL). 

We will use the following formula from IZaichik et al\ ( 20061 ) to connect the velocity 
scalings in different subranges, 



Su{l)^2u''- 



1 — exp 



4/3 



l^ + {2/CfL^ 



1/6 



(2.30) 



(15C)3/4,y 

The transverse structure function Snn then follows from eq. (2.26). 

The timescale T/j as a function of the separation in eq. ( 2.23) has d i fferen t scalings in 
the three subranges above as well. The theoretical model bv lLundgrenl (|l98ll ) found that 
in the viscous range, 

Tr{1) = VhTri, for I ^ (2.31) 

which was later confirmed by numerical simulations ( Girimaii fc Popdll990r) . 
In the inertial range, the similarity argument suggests that, 

Tr{1)^C21-^'^P'\ for r, <, I <, Li. (2.32) 

Following I Zaichik et ali ( 20061 ) . we will take the coefficient C2 ~ 0.3 in our calculations. 
Eqs. (2.31) and (2.32) connect at ~ (\/5/C2)^/2??, which is ~ 20r; for C2 = 0.3. This is a 
little larger than the corresponding transition scale (13??) for Su. 

For I :s> Li, the flow velocities across I are independent, Tr is thus expected to be the 
correlation timescale of the velocity along the trajectory of each p article, which is approx- 
imate ly given by the Lagrangian correlation timescale, Tl (see IZaichik fc Alipchenkov 
200i), 

Tr{1)^Tl, for l^Li. (2.33) 

This connects to the inertial-range scaling at about €:^^^{Tl/C2)'^^'^, which is ~ i using 
C2 = 0.3 and Tl ~ 0.3M'Ve. 

Similar to the case for Su, a formula is used to connect Tr in different subranges, 



^^(0 = Tl 



1 — exp 



V5j 



3/2 



-2/3 



l^ + 



Tl{e)yCl 



1/6 



(2.34) 



which is again adopted from lZaichik et all ( 2006[) . 

By analogy to the definition of the Taylor micro timescale, tt, we estimate ttr{1) 
by T^j^ ~ 2{5u{lY) / {6a{l)'^) where 5u{l) and 5a{l) are the velocity and acceleration 
difference across a distance I. This formula is simply a generalization of tt for the one- 
particle Lagrangian correlation to that for two-particle Lagrangian structure function. 
The correlation length scale of the acceleration field is expected to be short, probably 
~ 77. In that case, for I in the inertial range {5a{lY) ~ 2a? with a the rms acceleration 
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and thus ttr oc Su{l) - Z^/^. From Tr{1) - P^^, we have ttr{1) oc TrUY/^. In our 
calculations, we will use ttr — tt{Tr/Tl)^''^. Note that t his formula gives ttr — tt 
for I ^ L, as expected. It gives a constant ttr for / — >■ 0. IZaichik et al\ (|2003i 120061) 
assumed that ttr{1) = {tt /Tl)Tji{1) = zTr{1) without providing a physical motivation. 
Our calculations find that the two different assumptions for ttr do not give significant 
difi^erence in the predicted relative speed. 



2.3.2. Particle Pair Dispersion 

We now consider the rms distance, p(t), of two particles as a function of time r, which 
is needed to evaluate R (eq. (2.24)). As mentioned earlier, the specific question we ask 
here is how particles separate from each other backward in time, given their separation, 
r, at time zero (we will also refer to r as the initial distance for the backward dispersion 
from the viewpoint of the reversed time direction, although with normal time direction 
it is the final distance of the two particles in question) 

The study of turbuleiit disp e rsion of iiiertial p a rticles star t ed only recently (jBec et 



2007t iFouxon fc Horvai 1 120081: iBec et aZ.1 12009 



pa rticles star t ea only recently l|-Dec et at 
a|) , iBec et all (J2009al) gave a detailed re- 



port of simulation results for the separation behavior forward in time. They found that 
there are two temporal regimes with different separation behaviors: a transient regime 
representing the relaxation of the particle velocity toward the flow velocity, and a later 
regime where the particle pairs separate in a similar way as tracer particles. The transient 
regime lasts for about a friction timescale. For 5*^ ^ 3, a ballistic separation is found in 
the transient regime and the separation speed is equal to the initial velocity difference. 
The ballistic separation is due to the particles' memory for a period of ^ Tp. In the later 
phase, the particle separ ation is found to fo llow the Richardson- Obukhov separation law. 
Although the study bv lBec et all (|2009al) is for dispersion forward in time, their results 
provide a very useful guideline for us, because we are are not aware of any investigations 
for the b ackward d i spersio n of inertial particles. A separation behavior similar to that 
found bv lBec et all (|2009a|) will be used in our calculations (§3.1.3). We will show that 
a combination of an earlier ballistic phase and a later tracer-lik e phase gives quite good 
fit to the numerical simulation results for the relative speed bv lWang et all (|2000( ). For 
St <, 1, the particle separation appea rs to increase slow er than linearly with time in the 
transient phase according to Fig. 5 in lBec et aL ( 2009a ). Howev er, no function fit to the 



separation behavior in this regime is given bv lBec et all ()2009a|) . For simplicity, we will 
assume that the separation is ballistic for the early phase of all particles. This assumption 
gives rise to uncertainty in our prediction for the relative speed for particles with St <, 1. 

In order to understand the results of our model with the two-phase separation, we need 
to know the effect of each phase, and thus we first consider two simplified cases assum- 
ing a complete ballistic behavior (§3.1.1) and a complete tracer-like behavior (§3.1.2), 
respectively. The simplified cases give a very useful illustration for the effect of particle 
separation on the predicted relative speed. For example, for both ballistic and tracer- 
like separation behaviors, the relative speed as a function of the Stokes number can be 
explained from an approximate analysis of the integral equation for Vij . In particular, 
the analysis gives physical insights on the scaling of the relative speed with the Stokes 
number in the inertial range for the monodisperse case. Once the two simplified cases are 
understood, it is straightforward to interpret the prediction of our model with the more 
realistic separation behavior by a combination of two phases (§3.1.3). 

For ballistic motions, the particle distance goes linearly with time. Given the particle 
distance, r, at time zero, the separation as a function of r is. 



pHr) 



(2.35) 
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where [uP') = ((d^^^ — v'^'^^Y) — Spii{r) is the 3D relative velocity variance of the two 
particles at time zero. The separation speed at any time is taken to be the same as that 
at time zero. Recall that in our notation p(t) is the rms of p{t). The relative velocity 
variance, (ui^), in eq. (2.35) is unknown and is directly related to the radial relative speed 
under pursuit. We will build an implicit equation for {w"^) and (w^) in §3, which is then 
solved self-consistently. 

The exponential cutoffs in the integrand of eq. (2.14) imply that the primary contri- 
bution to Vij is from — t^i ^ t ^ and — Tn? ;$ r' ^ 0. This means that if the ballistic 
behavior lasts for about r^ ( Bee et aLll2009al) . then using eq. (2.35) at all times when 
integrating eq. (2.14) may give an acceptable order-of-magnitude estimate for T^ij, even 
though the separation is not ballistic at later time. 

Williams fc Grand (|l983l ) (and also lKruis fc Kusterslll997 ) take the particle motions to 



be ballistic in their calculations for the limit with large particles {Tp 3> T^). Assuming the 
velocity correlation between two particles can be neglected for the purpose of estimating 

1 /9 

the particle separation, they set the linear separation rate to ((u'^^')^ + (w''^^^)^) . We 
note that this assumed separation rate may be a good approximation only for very large 
particles. In the case of small to intermediate particles, the ve locities of nearby particle s 
are correlated and the separation rate is smaller than given bv lWilliams fc Grand ( 1983 ). 



Using their separation speed for those particles would overestimate the relative speed 
because both Sij and $ in eq. (2.23) increase with the particle distance, as can be seen 
from eqs. (2.30) and (2.34). 

In the second simplified c ase, we will consider the separation behavior sim ilar to that 
of tracer particle pairs (^see iFalkovich et al. II2OOII ISalazar fc Golling|2Q09l for detailed 



reviews on the pair dispersion of tracers). We are particularly interested in the effect 
of Richardson-Obukhov separation (the R-0 separation hereafter) phase found in the 
forward dispersion of inertial particle pairs. The R-0 separation law is written as 

p2(r)cxge-|T|3 (2.36) 

where the dimensionless coefficient g is known as the Richardson constant. iBec et all 
( 2009a[) did not give best- fit values for g in th e late-phase separa tion. Apparently g needs 



to be adjusted to fit the simulation results ()Bec et all (|2009al )). and it probably has a 



Stokes number dependence. Since the value of g for inertial particles is unknown, we will 
first use g measured for tracer pairs as a refe r ence. For tracers, theoretical models and 



direct numerical simulations bv lSawford et ali (|2Q05I ) show that the dispersion backward 



in time is sign i ficant ly faster than the forward dispersion. Experimental measurements 
bv lBerg et ali (|2006l ) found that g = 0.55 for the forward dispersion, and g = 1.15 for 



the dispersion backward in time. We will take g = 1.15 as a reference value since it is the 
backward dispersion that is relevant in our problem. In §3.1.3, we show that a two-phase 
separation with g '^ 1 in our model gives good fit to the numerical results for the relative 
speed. 

When the particle separation becomes larger than the integral scale, the flow velocities 
"seen" by the particle pair are uncorrelated and the separation is expected to be diffusive. 
Thus we will switch from the R-0 law to the diffusive separation when the particle 
distance exceeds L, i.e., we set /ci^(t) ~ 2Z?|r| for p ^ L where the coefficient D is given 
hy D = Qu'^Tl ■ Note that D here is for the 3D diffusion, and it is twice larger than D for 
1-particle diffusion (from the Taylor theorem) because the rms relative velocity of two 
faraway particles is 2u'^. We find that the exact separation behavior in the range p^ L 
in the diffusive regime is not important for the integration of eq. (2.14) because Sij and 
Tfl in eq. (2.23) becomes constant in this range of the particle distance. 
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The smallest initial separation considered bv iBec et al\ ( 2009al ) is about Irj. For an 



initial separation much below 77, there probably exists an initial exponential separation 
phase (similar to that of tracer pairs at small separations). The Ly apunov ex p onent s 
for particles with St f^ 2 have been computed from simulation by iBec et all ( 20061 ) . 



Although the exponential regime woul d also exis t for la rger particles as implied by the 



chaoticity of the dynamics (see, e.g., iBec et al\ (|2007r ) for a theoretical model which 



predicts that the Lyapunov exponent decreases as St~'^^^ for large St), its relevance at 
finite (but below 77) scale separation is questionable at least in the large St limit. Here, 
we will not consider the exponential phase, since it is unknown how long it lasts and 
how it connects with the lat er phases. In most o f our calculations we will give results for 
r ~ 77. In the simulations bv lWang et all ( 20001 ) and Zhou et all (|2001 ) that we will use 



to test our model, the relative speed is measured at a distance of 77, thus it is sufhcient 
to use Bee et al.^s result as a guideline in the comparison with those simulation results. 
It is straightforward to incorporate an exponential phase into our model, and once the 
detailed separation behavior of all particles at distances well below 77 is known, our model 
can predict the relative speed at any separation. 

Based on our physical picture, a careful consideration of the particle pair dispersion is 
necessary for an accurate estimate of the relative velocity. It is likely that the physics of 
turbulent dispersion of particle pairs is implicitly incorporated into the quite successful 
model by Zaichik and collaborators. In fact, their equation for the joint pdf of the particle 
separation and the relative velocity can be regarded as one for the pdf of the particle 
distance when integrated over the relative velocity (phase) space. For a comparison of 
our results with Zaichik et al. 's model, it is useful to see how tracer pairs separate in their 
framework, i.e., the prediction for particle separation from their fornrulation in the limit 
Tp — )■ 0. In this limit, the joint pdf equation can be reduced to an equation for the pdf of 
th e particle dis t ance. The pdf equation turns out to be in the same form as that suggested 
by RichardsonI (1 19261) . With structure functions and timescales given in §2.3 (same as in 






[Zaichik et al. (2OO67), the equation suggests that the Richardson constant, g, is about 3, 
significantly larger than measured from experiments and numerical simulations. It is not 
clear whether and how the quasi-normal assumption made by Zaichik et al. to close the 
moment equations of the joint pdf equation may affect the particle dispersion, or what 
it physically corresponds to regarding the separation. 

2.3.3. Average over the Direction of R 

We calculate the average of Sij over the direction of R in eq. (2.23). For R at any given 
time, we define a separation difference AR = R—r, the change of the separation from r at 
time zero. The direction of the separation difference is expected to be completely random 
if the flow velocity is statistically isotropic. This means that (Ai?iAi?j)ang = ^^^ij- 
We set I in eq. (2.25) to be ii = r + AR and take the average over the direction of AR. 
A rigorous derivation of this average needs to consider the dependence of Sij {R) on AR 
through the length R = \r + AR\ (e.g., in Su and 5„„) and that through the tensor 
RiRj simultaneously. However, the derivation is very complicated and cannot be done 
analytically. For simplicity, we neglect the dependence through R in the averaging process 
and only consider the average of RiRj over the AR direction. With this approximation, 
we find. 



{Sij{R))iine — 



2 ^2 \ /I r^ \ 

3 + 3^)^„„(i?)+(3-3^)^n(i?) 



% + (5„(i?)-5„„(i?))^ 



(2.37) 
where we have used (Ai?irj)ang — and {AR)^ — R^^ — r^. 
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The generalized shear term, I?y, then follows from eqs. (2.14), (2.23) and (2.37), 



_ , , f° dr f" dr' 
^vir)= — — 

> ^pl J~oc ^p2 



2 J^\ /I -2 



,j,.jSn4R)+[-,-^]SuiR) 



Stj 



+ {SuiR)-Snn{R)f-^\xHr-T';TTB{R),TR{R)) (2.38) 

f t\ (t' 

X exp — exp — 

\TplJ \Tp2 

We will numerically integrate this equation in §3 using the structure functions Su and 
Snn and the timescales ttr and Tr given in §2.3.1, and particle pair separation laws for 
p and R given in §2.3.2. 

In eq. (2.38), the structure functions, the timescale Tr, and hence $ (eq. (2.16)) in- 
crease with R, which increases with \t\ and \t'\. Together with the exponential cutoffs, 
this suggests that the integrand in eq. (2.38) peaks at r ~ — r^i and r' ~ —Tp2- Therefore, 
the main contribution to the integral is from r ^ — Tpi and r' ^ —Tp2 (if Tpi and Tp2 are 
not very different) and an important factor to determine the value of the integral is the 
distance, i?, at r ~ —Tpi and r' ^ — Tp2, which we will refer to as the primary distance. 
If Tpi and Tp2 are both very small so that the primary distance is close to r, Vij would 
approach Sij (r) as expected for the S-T limit. If the primary distance is much larger than 
the particle distance at time zero, r , the r^ /R^ terms in eq. (2.38) can be neglected and 
Vij ex 5ij (meaning that Vu = X>„„ = Vii/i). 

If Tpi and Tp2 are in the large limit and the primary distance is much larger than the 
length scale L, then typically the structure tensor in the integral in eq. (2.38) is ~ 2u'^dij 
(see eq. (2.29)), and the <& term would be the same as eq. (2.16) because Tr — Tl (eq. 
(2.33)) for R> L. This means that the integrand of eq. (2.38) would be the same as that 
in the 2nd term on the r.h.s. of eq. (2.13) (with B^ij given by the assumptions in §2.2) 
for the range of r and r' with R{t,t') <; L. On the other hand, at smaller values of \t\ 
and \t'\ (with R < L), the integrand for V^j is smaller than that in the latter. Thus 1?^^ 
given by (2.38) does not exactly cancel out the 2nd term on the r.h.s. of eq. (2.13). We 
find both numerically and analytically that the sum of the two terms decreases faster 
than 1/Tp in the limit Tp -> oo (it goes like Tp ' for the assumption of the ballistic 
separation and like ~ r"^ for the tracer- like separation behavior). Therefore, in the limit 
of large friction timescales our model satisfies the constraint discussed at the end of §2.1, 
namely, a good model needs to give a particle velocity correlation that decreases with Tp 
faster than the other terms. 

As discussed earlier, we used the rms distance, p and R, in our estimate of the trajec- 
tory structure tensor, S'tjj, while a rigorous derivation needs to consider the probability 
distribution function of the separation and take the average of S'tjj over this distribu- 
tion. Here we give a justification for this approximatioryj. First, if the separation is in 
the viscous range, R <,rj, Sij has a quadratic dependence on the separation (eq. (2.27)) 
and $ is independent of R (eq. (2.31)). This suggests that the approximation is exact 
since the average of the quadratic dependence over the separation distribution is ex- 
actly the square of the rms distance. If R is in the inertial range, then Sij ex R^'^, and, 
roughly speaking, the $ term provides another factor of R"^/^ (from the timescale Tr{R), 
eq. (2.32), over which the flow velocity difference "seen" by the particles is correlated). 

f In the justification argument, R is taken to be a stochastic variable, the length of the 
stochastic vector R, although at all other places in the paper R refers to the rms length of R. 
The rms distance will be written explicitly as (7?^)^" in the argument here. This change of 
notation is only for the present paragraph. 
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Therefore, using a single rms distance to replace the distance pdf would overestimate 
Stz] (and hence V^j) by a factor of ~ {R^)"^/^ / {R^'^). If the R pdf is 3D Gaussian as 
in the relative di ffusion model by Batchel or (1952), this factor is only 1.07. Numerical 
simulations (e.g., Bof fetta fc Sokolovll2002r) show that the separation pdf of tracer pairs 
is highly non-Gaussian with a very broa d tail. The pdf is found to be well fit by the 
solution of the pdf equation proposed bv iRichardsonI ()l926t ). With this broader pdf, we 
find the factor is larger, ~ 1.20. Because 5^ and Tr are independent of Rioi R^ i, our 
approximation is also expected to be exact for R in that range. In conclusion, replacing 
the R pdf by a delta function at the rms distance is quite well justified. It may overesti- 
mate Snj and Vij by ^ 20%, or the relative speed by ~ 10%, if the primary distance is 
in the inertial range (which is the case for Tp in the inertial range). 



3. Results 

To calculate the radial relative velocity, we need the longitudinal particle structure 
function, Spu, which in turn requires An and Vu. From eq. (2.19), it is clear that An = 
Ann and 



All 



(r!2-^l)^ r!if]2 + (^l+^2)^ 



n, + m 



(f7i + f72) [ni + ni 



(3.1) 



We obtain Vu using the relation Vu = Vijr.irj/r'^. From eq. (2.38) for Vij, we have 



Vii{r) 



° dr f" dr' 



oo ''pi J — CO ''p2 
X 



1 2r2 
3^3^ 



Sii{R) 



2 2r^\ 
3^3^/ 

T 



*(r - r'; ttr{R),Tr{R)) exp ( — j exp ( — 



i>nn\R) 

' t' 



(3.2) 



We will need the 3D relative velocity under the assumption of ballistic particle separa- 
tion. In that case, we take the 3D relative velocity variance, (w^), to be S^a = An -I- Vn 
(see §2.2). Using contractions of eqs. (2.19) and (2.38), we have, 

{w') = 3Aii+[ —[ ^{Sii{R)+2Snn{R))Hr-r';TTR,Tn)e^p(—)e^p(— 

J-ooTpl J-oo '^p2 K'^plJ \Tp2 

(3.3) 
where An is given by eq. (3.1). 



3.1. The Monodisperse Case 

As discussed earlier, for identical particles with Tpi = Tp2 — Tp, Aij vanishes and only 2?^^ 
contributes to the particle velocity structure tensor, i.e., Spij = Vij. Thus the longitudinal 
particle velocity structure function is given by. 



Spiiir) 



° dr r° dr' 



1 



2r' 



Sii{R) 



3 ■ 3i?2; "^ ^ ' V3 3i?V 

r 



2 2r^\ 



Snn{R) 



;$(r - r'; TTRiR), Tr{R)) exp - exp 



•p 



(3.4) 



To solve this equation, one needs the particle dispersion laws for p to calculate the 
distance R by eq. (2.24). We start with the two simplifi ed cases in order to study the 
effect of each phase in the two-phase separation found by iBec et all |2009a[ ) . 
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Figure 1. The radial relative velocity, (|wr|), as function of the Stokes number, St, for ballistic 
particle separation. The left panel shows results for different Rex. The distance r is fixed at 77. 
A St^" scaling is found between the S-T limit and the large St limit. The Stokes number at 
the transition from St^'"^ to St~^''^ scalings increases linearly with Re\. The right Panel shows 
the dependence on r with Re\ set to 300. 



3.1.1. Ballistic Separation Behavior 

We first consider the effect of the ballistic separation phase by using eq. (2.35) for the 
particle distance at all times. Eq. (3.3) gives the variance of the 3D relative velocity, 

(w^ 



), needed in eq. (2.35). For identical particles, eq. (3.3) becomes. 



/•o r]r r° dr' 

/ ^/ '^{Su{R) + 2Snn{R))Hl 

J - 00 ''p J -00 ''p 



T';TTR,TR)exp — 



exp [ — 

p 



(3.5) 

Because the distance R in the integral on the r.h.s. depends on (w^), eq. (3.5) is implicit 
for (w^). We numerically solve eq. (3.5) by an iterative method. After obtaining (w^), 
we use it to calculate Spu from eq. (3.4). 

The results for the radial relative velocity, (jwrl), as a function of the Stokes number, St, 
are shown in Fig. 1. We obtained {\wr\) from the conversion {\wr\) = y/2Spii/TT assuming 
a Gaussian distribution for Wr- We simply follow this convention here although we realize 
that the relative velocity distr ibution is broader than Gau ss ian for small partic les in high 
Reynolds- i iumber flows (e.g., ISundaram fc Collina (jlQQTl ): IWang et all (|2000l ). see also 
Bee et al\ ( 2009b[) for the intermittency in inertial particle structures). We will give the 
relative speed across fixed distances as a function of Stokes numbers (for convenience in 
the comparison with simulations in §3.1.3), although in coagulation models one needs to 
use the relative speed across the particle size as a function of the size. The latter can be 
easily calculated from our model with the Stokes number as a function of the particle 
size in a specific application. 

In the left panel of Fig. 1, we plot (|wr|) for particles at a distance r = ry as a function 
of the Stokes number and the Taylor Reynolds number. The solution reproduces the S-T 
limit (with the relative speed independent of St) and the St~^/^ scaling in the limit of 
large friction time. As in the model by Zaichik et al. , the relative velocity is found to scale 
as St^/"^ for intermediate Stokes numbers. This scaling corresponds to the inertial-range 
scaling of the turbulent flow and the St^^"^ scaling range will be referred to as the inertial 
range. 

The 5i^/^ scaling can be understood as follows. As discussed earlier, the main contribu- 
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tion to the integrals in eqs. (3.4) and (3.5) is from t', t ^ ~Tp, and the particle distance, 
R, at T, r' ~ — Tp, called the primary distance in §2.3.3, is important in determining the 
relative velocity. We will denote the primary distance as Rp. The primary distance as a 
function of Tp is evaluated by Rp = R{—Tp, —Tp) using cq. (2.22). For intermediate Stokes 
numbers, Rp is much larger than r, we thus have S-pu — (w^)/3 from eqs. (3.4) and (3.5). 
Thus the scaling behavior of {wf) with Tp is the same as that of {uP). The latter can be 
obtained by analyzing cq. (3.5). 

The $ term in eq. (3.5) represents the persistency of a structure of size R. It is ap- 
proximately given by exp(— |r — t'\/Tb,{R))- The effect of this factor depends on how 
Tr(Rp) at the primary distance, Rp, compares to Tp. If Tfi{Rp) is larger than r^, $ 
would be essentially unity for t and r' values that significantly contribute to the integral 
in eq. (3.5). One the other hand, if Th at Rp is smaller than Tp, a factor of Tr{Rp)/tp 
needs to be accounted because the $ factor suggests that, for a given t, only r' that 
satisfies |t' — t| <, Tn (instead of the range — Tp ^ r' ^ 0) contributes significantly to 
the integral. We find the latter is the case from our numerical solution. Thus, consider- 
ing the main contribution to the integral is from the integrand at R ^ Rp, we expect 
{w"^) oc [Sii{Rp) + 2Snn{Rp)\TR{Rp) / Tp. For balhstic separation, the primary distance 
Rp ~ {w'^)^pTp. And for R in the inertial range, we have Sii{R), 5„„(i?) ex R^^^ and 
Tr{R) cx i?2/3 (eqs. (2.28), (2.32)). These scalings give (w^) ex {w'^)'^/^Ty^ , which resuhs 
in (ui^) cx Tp and hence the St^^ scaling for the relative velocity. 

As the Stokes number increases, the relative speed reaches a peak, and transitions to 
the St~^P scaling. This change occurs when Rp ~ L and corresponds to the switch of the 
scaling behaviors of Su , S'„„ , and Tr from the inertial range to the outer scales. As pointed 
out earlier, if the primary distance, Rp, is much larger than L, the particle velocities 
are uncorrelated. A St~^P scaling is expected from the particle velocity variances, eq. 
(2.21). This scaling can also be obtained using an analysis similar to that for the Si^/"^ 
scaling given above. If Rp ^ L, the structure functions and the timescale Tr that give the 
primary contribution to the integrals in eqs. (3.4) and (3.5) are constant, i.e., Su = Snn = 
2u'2 and Tr ^ Tl (eqs. (2.29) and (2.33)). Therefore, eq. (3.5) gives {w^) cx u'^Tl/tp, 
and hence the St~^P scaling. Again the factor Ti^JTp comes from the $ term. 

The Stokes number, St^, or the friction time, Tpm, at which the relative speed reaches 
the maximum can be approximately obtained by setting Rp ~ L. Using Rp ~ {w'^Y''^Tp 
for the ballistic separation and {w^y^ ~ l.OSt^^Uri in the inertial range from our 
numerical solution, we find that Stm — -Rca/vIS, or Tpm — Te, where Tg is defined as 
u'^/e. The order-of-magnitude estimate for Stm turns out to be in good quantitative 
agreement with the numerical solution in the left panel of Fig. 1, which also confirms the 
linear increase of Stm with the -Re^. 

The right panel of Fig. 1 shows the dependence of the relative velocity on r, the particle 
distance at time zero. For r <, rj oi interest here, the relative velocity depends on r only in 
the S-T limit where it increases linearly with r (see the shear term in eq. (1.1)). For larger 
particles, it becomes independent of r. This is because in that case the contribution is 
mainly from the particle memory of the fiow velocity difference when the particle distance 
was much larger than r. 

We have used the same distance, R, for the timescale Tr in the $ term as that for 
the structure functions. As discussed in §2.2, for Tr it might be a better choice to use 
the particle distance, Rm, at the earlier one (min(r, r')) of the two times r and r'. With 
this choice, Tr is larger because the separation is larger at earlier time. The predicted 
relative velocity is also larger, and we find an increase by ~ 30% in the inertial range. 
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Figure 2. The radial relative velocity, {|wr|), at r — r] as a function of the Stokes number, 
St, assuming the R-O separation followed by a diffusive phase. The left panel shows results for 
different Rex- As in the case of ballistic separation, a St^'"^ scaling is found in the inertial range. 
The Stokes number where the curve peaks also increases linearly with Re\. The right panel 
gives the dependence of the relative velocity on g. In the inertial range it increases with g as 
«i/3 



3.1.2. Tracer-like Separation Behavior 

In this subsection, we examine the effect of the tracer-like separation phase. Our main 
purpose here is to study how the relative velocity scales with the Stokes number from 
the R-0 law. 

Here we will consider an initial distance of r = 77, and start the separation with 



p{r) 



,,2 



5e|Tp. We switch to the diffusive regime when the separation exceeds L 



(§2.3), and the connection is chosen such that p(t)^ = L^ + 2D\t — Td\ for |r| > |rd|, 
where Td is the time when p reaches L. In Fig. 2, we show the results for the relative 
velocity as a function of St and Re\ for particles a.t r = rj. The Richardson constant is set 
to be g = 1.15. The solution correctly reproduces the limits at small and large friction 
timescales as expected, and interestingly, we also find a St^^^ scaling for the relative 
velocity in the intermediate range of St. 

To explain the St^'-^ scaling, we analyze again Spu by considering the integrand in eq. 



(3.4) at the primary distance, Rp 



For Rp » r. 



eq. (3.4) can be approximately written 



as Spu ~ [Sii{Rp) + 2SnniRp)]TRiRp)/i3Tp). The factor Tr{R(tp))/tp (which is smaller 
than 1 for any g <^ 10) is based on the same reasoning as in §3.1.1 for the ballistic 
case. With the R-0 separation law, the primary distance is given by Rp ~ (geTpY^^. 
Therefore, using the inertial-range scalings of 5*;;, 5'„„, and Tr, we have Spu ex g'^^^Tp. 
This explains the S'i^/^ scaling for the radial relative velocity in the inertial range, and 
also predicts that the relative velocity increases with the Richardson constant as g^^^. 
The g^^^ dependence is confirmed by the right panel of Fig. 2. The increase of the relative 
velocity with g shows that faster particle separation gives larger relative velocity. 

As in the ballistic case, the friction time, Tpm, at which the relative velocity peaks, 
is again obtained by setting the primary distance Rp to L. Using the R-0 separation 
law for Rp, we find that Tp„i ~ g~^^^Tf,, with Tg = u'^/e. In units of the Kolmogorov 
timescale, we have St„i ~ g~^^'^Rex/VT5. This is consistent with the results in the left 
panel of Fig. 2. 

We have finished the study for the effects of each separation phase in the two-phase 
separation found bv lBec et all (|2009a[) . We found that, in the inertial range, the relative 
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speed from using the R-0 law has the same scahng behavior as from the balhstic sep- 
aration. Thus the same St^^"^ scahng is also expected for a combination of ballistic and 
R-0 separation behaviors. For g ^ 1, the predicted value for the relative speed in the 
inertial range from the R-0 separation is quantitatively very close to that in the ballistic 
case, and so is the peak Stokes number Stm- Therefore for the combined behavior the 
prediction for the inertial range and for the switch to large Stokes numbers would be 
similar to the two simplified cases (see Fig. 4). However, in the transition region from 
the S-T limit to inertial range, the predicted relative speed is quite different for the two 
separation behaviors. 

3.1.3. Combined Separation Behavior and Comparison with Simulation Results 

We use the simulation results bv lWang et al\ ( 200(1 1 to test our model. The Reynolds 



numbers i n their simu l ations are quite low with Re\ in the range from 45 to 75. The 
results bv IWang et all ( 20001 ) for three different Reynolds numbers are shown as data 



points in Fig. 3. The relative speed is for particles at a separation of rj. Due to the 
limited resolution of these simulations, no inertial-range scaling is seen. For Rex in the 
same range, our model with the two separation behaviors considered above does not show 
the St^^^ scaling either (see Rex = 100 curves in Figs. 1 and 2). 

The separati on behavior o f iner tial particle pairs at Reynolds numbers as those in the 
simulations bv IWang et al\ ( 20001 ) is unknown. We tried different separation behaviors 



and compared the predictions for the relative speed with their results. We find that a 
combin ation of an early ballistic phase and a later t racer-like phase, as found bv lBec et al. 



(|2009al ). can give a quite good fit to the results bv lWang et all (J2000l ) (while with a pure 



ballistic separation or a pure tracer- like separation no satisfactory fit is found) . The lines 
in Fig. 3 show the predicted relative velocities with such a two-phase separation, which 
agree quite well with the simulation data. The exact separation behavior used in Fig. (3) 
is as follows. The separation starts with a ballistic phase which lasts from time zero back 
to Tc. At Tc, it continuously connects to the R-0 separation law, and finally switches to 
the diffusive regime when the separation exceeds L. The connection between the ballistic 
phase (eq. (2.35)) and the R-0 phase (eq. 2.36) is chosen such that the particle distance 
at T < Tc is given by p{t)'^ = p(Tc)^-|-g|r — Tc|t^. This connection between the two phases 
is quite smooth. Furthermore, we set Tc = — 1.4rp and g ~ 1. The connection between 
the R-0 regime and the diffusive regime is the same as that for the tracer-like separation 
behavior discussed in §3.1.2. Due to the low Reynolds numbers here, there is only a 
very short period for the R-0 separation between the ballistic regime and the diffusive 
regime in the chosen separation behavior. How the ballistic phase exactly connects to 
the R-0 phase in the backward dispersion of inertial particles is unknown. We also tried 
other ways to connect the two phases, and found that, by adjusting Tc and g, some other 
connections can also give satisfactory fits. For example, if the two phases are connected 
by p{'t)'^ = p{tc)'^ + 9{t — '''c)'^\'''\ for r < Tc, fitting the simulation data gives Tc = — l.lrp 
and g — 1.5; and for a connection with /9(t)^ = pirc)"^ + qItI"^ — skcP in the R-0 phase, 
we find a good fit with Tc = l.Sxp and g = 0.6. 

Because Vij is an integral over the history of the particle distance, different separations 
as a function of time could lead to the same predicted relative speed. Therefore, the exact 
backward dispersion behavior cannot be determined by fitting the simulation results for 
the relative speed. In other worlds, a separation law that fits the data may not represent 
the ex act dispersion behavior for inertial particles in the simulated fiows bv lWang et al 



( 2000I ). and a verification would need a direct numerical study of the pair separation. 



However, the fact that the dis persion behavior used in Fig. 3 is generally consistent 
with the simulations results by iBec et all (|2009ah (see their Fig. 8 where the connection 
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St 



Figure 3. Comparison with the simulation results from IWang et al\ IJ2000I ) at Re\ = 45 
(squares), 58 (circles) and 75 (triangles). Lines are the predicted relative speed from our model 
with a two-phase separation behavior, a ballistic phase followed a tracer-like phase. The ballistic 
separation is assumed to connect with the R-O separation with g = 1 aX — 1.4rp. See text for 
details on this connection. 



between the ballistic phase and the R-0 phase occurs between 1-2 Tp ) suggests that the 
adopted behavior is at least qualitatively correct. We will adopt the separation behavior 
used in Fig. 3 for all the calculations in the rest of the paper. 

We find a significant deviation ( up to 25 %) between the model prediction and the 
simulation results for very large particles [St <; 20). It seems that the deviation at these 
large Stokes numbers could not be removed by a reasonable change in the separation 
behavior without causing discrepancy at smaller St. This deviation also occurs in Zaichik 
et al.^s model. An immediate suspect for this deviation is the assumption in our model 
(and in Zaichik et al.'s model) that the trajectories of all particles are not far away from 
those of the fluid elements This assumption is not well justified for very large particles. 
As discussed at the end of §2.2, the temporal correlation of the flow velocity (or the 
velocity difference) along the trajectories of large particles may be close to the Eulerian 
correlation, while we used the Lagrangian correlation timescale throughout the model. 
If the Eulerian correlation timescale were used for these particles, and if the Eulerian 
correlation timescale is larger than the Lagrangian timescale (see discussions in §2.2), the 
predicted relative speed for large particles would be larger, reducing the the difference 
between the model and the simulation results. Our approximation for the trajectory 
structure tensor in §2.2 could also contribute to the discrepancy. 

Using the separation behavior in Fig. 3 that well fits the simulation results, we carried 
out calculations for larger Reynolds numbers. Fig. 4 shows our prediction for Re\ = 300. 
For comparison, we included results assuming pure ballistic separation (the dashed thin 
line) and pure tracer-like separation (the solid thin line). For the latter we used 5=1 
for the R-0 separation. We see that, in the two-phase case, the relative velocity in the 
transition region from the S-T limit to the inertial range lies between the two cases with 
single separation behavior. In this transition region, the predicted relative speed relies 
on how the ballistic phase and the R-0 phase exactly connects. In the inertial range, we 
again have the St^^'^ scaling, which is expected because both the ballistic behavior and 
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Figure 4. The radial relative velocity, (j^yrl), as a function of the Stokes number, St, with the 
same two-phase separation used in Fig. 3. The dashed thin line is for ballistic separation only 
and the solid thin line is for tracer-like behavior only (with g = 1 for the R-0 separation). Also 
shown is the prediction by the Zaichik et a/.'s model. 



the R-0 separation law give that scahng. For g ~ 1, the predicted relative speed turns 
out to be very similar for all the three cases in the inertial range. The predicted speed 
in the inertial range is not significantly affected by the details of the connection between 
the ballistic and the R-0 phases. The Stokes number at the transition from the inertial 
range to the St^^^^ range is ~ Re\/\/T5 (r^ ~ Te), similar to both the pure ballistic 
case and the pure tracer-like case. 

Fig. 4 also shows the relative speed predicted by Zaichik et a/.'s model (the thin dotted 
line). We obtained the results of t heir m odel by numerically solving the set of differential 
equations given in lZaichik et all ( 20061 ) (i.e., their eqs. (51-53)). The transition from the 
S-T limit to the inertial range in Zaichik et al.'s model is quite steep. The relative speed 
in the inertial range in their model is 20% larger than the our model. This is probably 
because a faster particle separation is built into Zaichik et aVs model. As pointed our 
earlier, in the limit of Tp — >■ 0, Zaichik et aVs formulation implies a Richardson constant 
significantly larger than 1. 

3.2. The Bidisperse Case 

We first point out the importance of gravity in the bidisperse case. In the monodisperse 
case, the settling velocity by gravity is the same for all particles, and thus neglecting 
gravity in that case may give approximately good estimates for the relative speed. The 
situation is quite different in the bidisperse case where gravity may play a major role, es- 
pecially for large particles. The terminal velocity difference between two different particle 
could give substantial or even dominant contribution to their relative speed. Besides the 
direct contribution to the relative speed, differential settling can also have an indirect ef- 
fect by increasing the particle separation in the past. Since a larger separation backward 
in time would increases the contribution from the particle memory of the flow velocity 
difference, this indirect effect also tends to give a larger collision speed. 

We will neglect gravity in our model below. Clearly, this limits the application of 
the model only to situations where the relative speed caused by differential settling is 
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Figure 5. Comparison with the simulation results from IZhou et al\ l|200ll ) at Re\ = 45 
(squares), 58 (circles). Lines are the predicted relative speed from our model with a two- phase 
separation behavior, a ballistic separation followed by a tracer-like behavior. The ballistic sepa- 



ration is assumed to connect with the R-O separation with g ■ 
about the connection. 



1 at — 1.4rn. See text for details 



negligible in comparison to the prediction of our model. However, understanding the 
simpler case with turbulence alone is of theoretical importance because it serves as the 
first step to a physical and accurate model for particle collisions in realistic environments 
where both turbulence and gravity are present. The effect of gravity may be included 
in our framework by accounting for both the direct effect of differential settling and its 
indir ect effect through t he backward separation of the particles (e.g., in a similar way 



Avala et al\ ( 20081) . who, however, did not include the particle separation by the 



as m 
turbulent flow). 

In the bidisperse case, we include the contribution from the An term, eq. (3.1), to the 
relative speed. We use the same numerical method to solve Vu as in the monodisperse 
case. The particle separation behavior is chosen to be similar to the one that well fits 
the simulation results in the monodisperse case. Note that in the ballistic phase, the 
separation speed (i.e., the relative velocity variance, (w^)) to be used for the calculation 
oiDii has a contribution from the An term (eq. 3.3). It is not clear how long the ballistic 
phase lasts in the bidisperse case, since the friction timescales of the two particles are 
different. We simply assume that the duration of the ballistic separation is proportional 
to the average of the two friction timescales. Fig. 5 shows the radial relative velocity as 
a function of the Stokes number of particle (2) for a fixe d Stokes number, St^ = 1, of 
particle (1). The data points are simulation results from Zhou et al\ (20011 ) for r — rj 
at Re\ = 45 and 58. These data are from Fig. 15 of lZhou et al\ (J2001I ). but a different 
normalization is used here. The lines are the prediction of our model where we adopted 
the same connection used in Fig. 3 for the monodisperse case and set g — 1 and t^ = 
— 1.4 X (rpi -I- Tp2)/2. The agreement of the model prediction with the simulation results 
is quite good, except that it is a little bit broader around the dip. A possible reason is 
that the separation behavior for different particles is different from the one used here, 
which is based on the separation behavior of identical particle pairs. 

We give more details for the relative speed between different particles in Fig. 6, where 
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Figure 6. The radial relative velocity, (|wr|) as a function of St2 with fixed Sti. The thick lines 
give the relative speed at r = 77. Each curve has a dip at St2 = Sti. The three thin dotted lines 
show the results for r = 0.3r; with St — 1, 3, and 10. 



the Taylor Reynolds number is set to 300. The thick lines are for r — r]. A dip around 
5*^2 ^ Sti is found in every curve with Sti < 100. The existence of the dips is related to 
the fact that the contribution from Aij for particles of similar sizes is small. Physically, 
it means that the velocities of similar particles tend to have stronger correlation than 
particles with very different friction timescales. Around each dip, the contribution to the 
relative velocity is mainly from 2?n, while far from the dip it is dominated by Au. A 
comparison of the two terms shows that they give similar contributions when the Stokes 
number ratio is about 3-4. If the friction timescales of the two particles differ by a factor 
much larger than 4, using the generalized acceleration term alone may give a satisfactory 
result. The T>ii and Au terms in our mod el are close l y rela ted to the two terms in the 
equation for the velocity difference given in lBec et al\ (20051) (their equation (13)). Their 
discussion on the relative importance of those two terms provides physical insights to 
understand the dips in Fig. 6. The relative speeds at the dip centers correspond to equal 
Stokes numbers. Connecting these centers would give a curve identical to that for the 
monodisperse case with the same separation behavior and parameters. 

On the far left of the dip, the relative velocity approaches a constant. The constant 
corresponds to Au in the limit Q2 ~^ 0, which, from eq. (3.1), is given by „ +n'i\ 2/2 " ■ 
As Sti (or ill) increases, the dip moves to the right, and the relative velocity on the far 
left increases. It reaches and stays at the maximum (corresponding to Au = w'^) after 
ill becomes much larger than 1. The opposite occurs on the far right of the dip, i.e., in 



In this limit, Au approaches 






which has a maximum 



the limit St2 > 5ti. 

of m'^ at small Sti. With increasing Sti or ili, the relative velocity decreases on the far 
right of the dip, while it increases on the other side. 

The dip disappears for very large Sti (larger than ^ 100 for the case shown in the 
Fig. 6). This can be explained as follows. Physically, the dips are due to strong velocity 
correlation between particles of similar sizes. Thus no dip would exist if the friction time 
of particle (1) is such that its velocity is not significantly correlated with any particle of 
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similar size. This is the case for a particle with Tp larger than Tprm the friction timescale 
where the relative speed peaks in the monodisperse case (see §3.1). The velocity of such 
a particle is not correlated even with an identical particle. Therefore, no dip is expected 
around St2 — Sti if Tpi > Tpm — Tf,. Clearly, this argument suggests that the critical 
value of 5*^1 where the dip starts to disappear is the same as the Stokes number where 
the relative speed peaks in the monodisperse case. 

The dependence of the relative velocity on the distance r for the bidisperse case is also 
illustrated in Fig. 6, where the dotted thin lines show results for r — 0.3r;. The r depen- 
dence only comes from !?«, because Au does not depend on r. Thus the dependence may 
exist only around dips where !?« gives a significant contribution. As in the monodisperse 
case, I?n is independent of r (again for r < -q) ii the Stokes numbers of both particles are 
much larger than 1. This suggests that, in the bidisperse case, the relative velocity is a 
function of r only when 5^2 — Sti ^ 1. In Fig. 6 we see that for Sti — 1 the depth of 
the dip increases with decreasing r, corresponding to the decrease of the relative velocity 
with r in eq. (1.1) for small identical particles. The r-dependence is already weak for 
St\ = 3, and the relative velocity becomes completely independent of r for Sti = 10. 

In summary, we find that in the bidisperse case the dominant contribution to the 
relative velocity is from Dn if the ratio of the two Stokes numbers is not larger than 3-4, 
or from Au if that ratio is larger. For a fixed Stokes number of one particle, 5ii, the 
relative velocity as a function of the Stokes number of the other particle, 5^2, shows a dip 
at 5^2 ^ Sti. The dip corresponds to a stronger velocity correlation between particles of 
similar sizes than between very different particles. The ex istence of a dip at equal Stokes 
numb ers has consequences for the collision kernel (e.g., IBec et al\ l2005(: Zaichik et al\ 
200i). 



4. Comparison with Other Models 

Our model has already been compared with some previous models earlier in the paper. 
Here we discuss more models and present the comparisons more systematically. 

4.1. Volk et al. (1980) 

In the astrophysical literature, the relative velocities of dust grains, e.g., in the context 
of the coagulation growth of dust grains in protoplan etary disks, has almost always been 



01 tne coagulation growth ot dust grams m protoplan etary disks, has almost always been 
estimated b ased on the model bv IVolk et al\ (1980|), and a later version of that model 



bv iMarki cwicz et al\ ([1991). Here we discuss this model for the particle relative velocity 
and compare it with the physical picture presented in this paper. 

Volk et al (1980) derived both the 1-particle rms velocity and the relative speed be- 
tween two particles. Their model started by considering the effects of turbulent eddies of 
different sizes on a particle of a given friction timescale. They speculated that the effect 
of an eddy on the particle depends on the eddy size. The particle would basically move 
along with large eddies if the eddy turnover time is much larger than the particle friction 
timescale. On the other hand, the effect of eddies with turnover time much smaller than 
the friction time is argued to be like a "random kick" because the eddy would "die" 
within a friction timescale. Apparently "random kick" here only means that the particle 
does not "receive" the driving by these small eddies to a full extent. There is also another 
way that an eddy may behave like a random kick. Due to their inertia, particles have a 
different velocity from the flow. If the relative velocity between a particle and a turbulent 
eddy of a given size is such that the particle crosses and leaves the eddy within a friction 
timescale, then the particle would not receive a "full" kick from that eddy, and the effect 
of the eddy is a "random kick". Based on these considerations, Volk et al. defined a 
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critical eddy size in Fourier space, k* , such that the effect of eddies below this size is like 
random kicks. 

In the model by Volk et al. (1980) for the rms velocity of a single particle, the wavenum- 
ber k* appears to be important because eddies of size smaller than the scale I*, corre- 
sponding to k* , are expected to be less efficient at "driving" particle motions than larger 
ones because the particle does not have chance to "fully" receive the energy from these 
eddies. 

Volk et aZ.'s derivation for the relative speed between two particles did not consider the 
separation of two particles, which is essential in our physical picture. In their calculation 
for the relative speed, they continue to take k* as a crucial scale. They essentially assumed 
the velocities of two particles induced by eddies with wavenumber larger than k* are not 
correlated, and the contribution to the particle velocity correlation is only from larger 
eddies. This assumption is not justified. 

We argue that it is unlikely that k* is crucial for the velocity correlation of two particles 
(although k* may be an important scale for the 1-particle rms velocity, as discussed 
above). Whether the particle relative motions induced by eddies at a given scale are 
correlated or not is probably determined by how the particle separation compares to 
the eddy size. If the distance of two particles is smaller than the size of an eddy they 
encounter, the particle motions induced by this eddy would be correlated even if the eddy 
size is smaller than /* (the contribution to the two particles from a same eddy should 
be correlated). In this case, Volk et aZ.'s assumption would underestimate the correlation 
and overestimate the relative velocity. On the other hand, contributions to the velocities 
of the two particles from eddies of size smaller than the particle distance would be 
uncorrelated, because each particle receives a contribution from a different eddy of that 
size and motions in different eddies are likely to be independent. Therefore, contrary to 
the assumption by Volk et al, particle motions induced by eddies of size larger than I* 
are not always correlated. They are independent if the size of these eddies is smaller than 
the particle separation. The argument above suggests that it is the particle separation 
(instead of I* or k*) that determines the particle velocity correlation, and thus an explicit 
examination of the particle separation is required. 

The model by Volk et al. for the relative speed for two particles may be interpreted as 
one that implicitly assumes that the typical particle separation, R, is around the scale 
k*. However, this assumption cannot be physically justified since the definition of k* has 
nothing to do with the distance between two particles. Even if the value of k* turned out 
to be close the to typical distance between particles, it should probably be taken as a 
coincidence. Volk et a/.'s model may be improved by incorporating the particle distance 
as a function of time within their formulation. 



4.2. I Williams & Cran^ 1(198 A) 



Williams fc Grand (|1983I ) started from the derivation of the relative velocities for particles 
in two limits: Tp ^ Tl and Tp ^ T^, and then obtained a "universal" formula by 
interpolation. In the limit of Tp -^Tl, they assumed that the particle separation back in 
time can be neglected (§2.2), and found the relative velocity is given by eq. (2.20). As 
argued in §2.2, this result is not valid for similar particles especially for St ^ 1. In the 
other limit of Tp ^ Tl, Williams and Crane (1983) considered a linear particle separation, 
and chose the separation rate assuming the particle velocities are not correlated in this 

li mit. 

Williams fc Grand (|l983l ) found a "universal" formula that reproduces the results for 



two limits they considered. This formula is obtained from a mathematical interpolation, 
and thus does not incorporate the physics of the relative velocities between particles of 
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intermediate inertia. We find that the formula gives a St^^'^ scaling for identical particles 
with Tp ^ Tl, which is probably incorrect. 



4.3. [^ h98A) 



Yud (|1984[ ) derived a formula for the the relative velocity that consists of an acceleration 
term and a shear term. Neglecting the effect of the added mass term (which is negligible 
in gaseous flows) included in Yuu's calculations, the acceleration term is also given by eq. 
(2.20), i.e., exactly the same as Williams and Crane's result for the small particle limit. 
Yuu's shear term is much smaller than our generalized shear term, T^ij, for intermediate 
to large particles, because his calculation did not consider the particle separation back 
in time, and thus the shear term does not account for particles' memory of the larger 
flow velocity difference they "saw" at earlier times. As pointed out in §.2.2, this leads to 
a significant underestimate of the contribution from Vij. As a consequence, Yuu's model 
is not valid for similar particles with Tp'> Trj. 

4.4. \Krms & Kuste^ h99i) 



Kruis fc Kusterd (|l997l ) gave a generalization to the models of lWilliams fc Crana (jl983() 



and IYuuI (|1984I ). They first noticed that Williams and Crane's result for particles with 



Tp <C Tl and Yuu's acceleration term (eq. (2.20)) do not repro duce that in e g. (1) for 



the S-T limit. Replacing the temporal energy spectrum in Williams &: Crand (|1983l ) by 



one that incorporates the flow acceleration and corresponds to a temporal correlation 
function similar to our eq. (2.16), they were able to derive a formula for the Tp <^ Tl 
case that correctly reduces to the accelerati on term in the S-T limit . The formula is the 
same as our result for Aij in the limit z <^ l. lKruis fc Kusterj (|l997h also generalize d the 



Williams and Crane model to include the added mass effect, which was considered in lYuu 



( |l984 ). The effect is negligible in a gaseous flow, but may be important in liquid flows. 
Following Williams and Crane, a "universal" solution was obtained by interpolating the 
generalized results for t he two limits. Therefore the Kruis and Kusters model shares 
the same weakness as W illiams fc Cran e (1983): the physical importance of the particle 
separation for the relative velocity of similar particles with Tp <, Tl is not included in 
the model, and the interpolated results for that case are probably incorrect. 

4.5. The analytical model of Zaichik et al. (2003, 2006) 

In addition to their differential model discussed earlier in details, Zaichik and collab- 
orators also presented an analytical model. Assuming Gaussian statistics for both the 
flow and the particle velocities, the analytical model calculates the joint pdf of the ve- 
locities of two particles, P(t;^^\t> ^^^), from the two-point joint pdf of particle and flow 
velocities, P(t>(i).t)(2).M(i),tx(2)v |Zaichik et oil (2003. 2006) approximated the latter by 
P(t,(i) |u(i))P(«(2) |ii(2))P(u(i), m(2)) under the assumption that P(t;(i) |t;(2), «(!), ^(2)) = 
P(t!(i)|M(i)) and P(t;(2)|ti(i),M(2)) = P(t;(2)|n(2)), This assumption is valid in the limit 
Tp — > 0. The particle velocity is well approximated by the flow velocity at the same 
point in this limit so that P(t)(^^|t)(2) ^,(1) .j^(2)') ^ (5(^(1) — t>(2))^ Therefore it can be 
approximated by P(t>(^)|?;(2))^ which is also 5{v^^'^ —v'^'^^) in the limit. A similar argument 
applies to P(t;(2)|M(i)^M(2)) ^ P(t)(2)|M(2)). in the other limit with Tp > Tl, wc find the 
assumption for the conditional pdfs is also roughly valid. In this limit, due to the long 
memory, the particle velocity, v^^' is not strongly correlated with the local flow velocity 
m(i), nor with u^^) or m(2). Therefore P(u(i)|ti(2),M(i),ii(2)) = P(t;(i)|M(i)) could be a 
good approximation in the large particle limit. (The predicted relative velocity by this 
model could reproduce the two limits given in §1). However, for Tp in the inertial range, 
all the four velocities are partially correlated. The assumption that neglects the "direct" 
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correlation of the particle velocity at point (1) with the flow and particle velocities at 
point (2) would, to some degree, underestimate the particl e velocity correlation . 

From the 2-point joint pdf of particle and flow velocities. IZaichik et all ( 20031 ) derived 
the joint pdf P{v'^^\v'^^^) and an analytical formula for the relative velocity. In the 
derivation, the Lagrangian correlation function is needed to calculate the velocity variance 
of each particle and the flow-particle correlation at each point, which fix P{v^^''\u'^^'>) 
and P{v^^'>\u'^^'>) under the assumption of Gaussian statistics. With our eq. (2.16) for 
the Lagrangian correlation function, the model gives a linear scaling with St in the 
inertial range for the monodisperse case. This is in contrast to the St^^^ scaling in both 
their differential model and our model. Because the analytical model underestimates the 
velocity correlation, the predicted relative velocity for identical particles in the inertial 
range is much larger than in the latte r two models, as w ell as than the results from 
simulations with low Reynolds numbers ( Zaichik et aZ.1120031 ). Furthermore, the analytical 



model does not give a dip at 5*^2 — Sti in the bidisperse case. This is again because the 
velocity correlation between similar particles is not accurately evaluated. The model does 
not account for the fact that similar size particles tend to have stronger correlation than 
different size ones. 



4.6. iDerivicli JEOOd ) 



The model of lPerivichl (20061) starts with a similar approach as in the differential model 



of Zaichik et al. An equation for the joint pdf of both positions and velocities of the two 
particles is derived assuming Gaussian statistics for the flow velocity field. The equation 
is equivalent to that for the joint pdf of the particle separation and the relative velocity 
in Zaichik et al.'s model. Apparently, Derivich (2006) considered the solution of the joint 
pdf equation only in the zero separation limit pj. Neglecting the spatial derivative terms 
in the joint pdf equation (i.e., assuming the spatial derivatives can be neglected in the 
zero separation limit), Derivich obtained a solution for the joint pdf (Gaussian) of the 
velocities of the two particles. The correlation of the particle velocities in the solution 
depends on the particle separation as a function of time (with initial separation set to 
be zero). In his calculations, the particle separation is taken to be a Gaussian variable. 
The time dependence of the separation variance is neglected, and apparently for the 
monodisperse case the variance is set to be a constant, corresponding to the particle 
distance at r, r' — —Tp in the ballistic separation behavior assumed in our model. 
We argue that this treatment with a constant particle separation variance is physically 
inadequate. 

4.7. Ayala et al. 2008 
Our fo rmulation is very similar to that in lAvala et al\ ( 20081 ). The model bv lAvala et al. 



(|2008l ) included the particle separation due to gravity for sedimenting droplets in turbu- 



lent flows, but neglected the particle separation by turbulent dispersion. In the absence 
of gravity, particles do not separate in their model, and the model enormously underesti- 
mates the relative velocity in the monodisperse case (where the separ ation plays a crucial 
role), as can be seen from their Fig. lib. As explicitly pointed out bv lAvala et all (|2008l ). 



their model was not designed for particles of similar sizes. In the presence of gravity, it is 
expected the accuracy of the model decreases with increasing turbulence intensity. This 
can be seen from their fig. 13 for the monodisperse case, where the predicted relative 
velocity agrees with the simulation results for the lower of the two turbulent intensities 

f Zaichik et al. considered the whole range of the particle separations when solving their joint 
pdf equation and thus were able to examine the spatial clustering of particles that Derevich's 
model does not address. 
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shown (e — 100 cm^/s'^), while it is significantly smaller than the simulation results for 
the case with the higher intensity (e — 400 cm^/s^). Clearly, with larger turbulent in- 
tensity, the turbulent dispersion is faster and neglecting it would result in less reliable 
predictions. We also note that Ayala et al. adopted a bi-exponential form (similar to the 
form of our eq. (2.16) for the Lagarangian temporal correlation function) for the spatial 
correlation function of the flow velocity. The form corresponds to a linear velocity dif- 
ference scaling and hence a —2 energy spectrum in the inertial range (for comparsion, 
see eq. (2.29) for the spatial correlation function adopted in our model), and is thus not 
consistent with the Kolmogorov spectrum observed in turbulent flows of high Reynolds 
numbers. 



5. Conclusions 

We have examined the relative velocity of inertial particles suspended in turbulent 
flows. A general formulation is established based on the calculation of the particle velocity 
structure function. Our general result for the particle structure function, eq. (2.12), has 
two terms, a generalized acceleration term, Aij (eq. (2.13)), and a generalized shear term, 
T>ij (eq. (2.14)). The generalized shear term, Pjj, corresponds to the contribution to the 
relative speed from particles' memory of the flow velocity difference in the past. We find 
that the backward-in-time dispersion of inertial particle pairs is needed to evaluate this 
term. The two terms reduce to the acceleration term and the shear term, respectively, 
in the S-T limit. Our formulation can thus be viewed as a generalization of Saffman and 
Turner's result for the limit of small particles to particles of any size. 

We have shown that our model with a separation behavior similar to that found b y 
recent simulations for the forward (in time) dispersion of inertial pairs (IBec et a/.ll2009al ). 



i.e., a ballistic separation followed by a trace r-like behavior, give s quite good fits to the 
relative spee d measured from s imulations by IWang et al\ (J200G| ) (for the monodisperse 
case) and bv lZhou et al\ ( 200ll ) (for the bidisperse case). 



For the monodisperse case, only Vij contributes to the relative velocity. At large 
Reynolds numbers, a St^^'^ scaling of the relative velocity in the inertial range is found for 
both the ballistic separation or the Richardson separation. Therefore, for the two-phase 
separation that well fits the simulation results, we have the same inertial-range scaling. 
This sca ling is consistent with that f r om the differenical model by Zaichik and collab- 
orators ( Zaichik fc Alipchenkovll2003t IZaichik et aZ.ll2003l 120061 ). Our model provides a 



clear physical picture for this scaling. 

Our calculations for the bidisperse case show that Aij dominates the contribution to 
the relative velocity between particles of very different sizes, while for similar particles 
the primary contribution is from "Dij. In the relative velocity vs. St2 curves with fixed 
5^1, dips are found around 5*^2 ^ •S'ii, indicating stronger velocity correlation for similar 
size particles than for different size ones. Away from the dips, the relative velocity is 
essentially given by the contribution from Aij . 

The main assumptions in our model are those for the trajectory correlation and trajec- 
tory structure tensors. The approximations for these tensors can be tested and improved 
by numerical simulations, and our work thus provides a motivation for direct studies of 
these correlations along the particle trajectories. A direct numerical study of the sepa- 
ration behavior of particles backward in time would also be of interest, because we have 
shown that it plays an important role in modeling the relative velocity between particles 
of similar sizes. With the help of future simulations, the assumptions in our model could 
be considerably refined. The model may also be extended to include gravity and other 
effects. The refined and extended model would provide a reliable prediction of the relative 
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velocity between inertial particles and can be applied to many practical studies, such as 
raindrop formation in atmospheric clouds and collisions of dust grains in astrophysical 
environments. 

LB acknowledges support from NASA grants 08-NAI5-0018 and NNX09AD10G. 
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